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ABSTRACT 


Our  goal  was  to  measure  the  distribution  of  radiation  S  as  a 
function  of  its  angular  coordinatesa,  g  and  its  wavelength  A.  We 
first  measured  S(ct),  i.e*  the  distribution  of  radiation  as  a  function 
of  one  angular  coordinate.  From  this  we  learned  how  critical  the 
various  parameters  of  our  instrument,  such  as  collimation  and  the 
gratings’  substrate  flatness,  were  to  its  successful  operation. 

These  investigations  led  to  the  design  of  new  interferometric  test 
instruments  which  include  an  autocollimator,  and  a  lateral  and  radial 
shearing  interferometer.  Finally  a  theoretical  study  of  an  instru- 
ment  to  measure  S(a,g,X)  was  accomplished. 
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FOREWARD 


In  1836  H.  F.  Talbot  discovered  that  images  of  a  grating  illumi¬ 
nated  with  a  point  source  of  light  were  formed  without  the  aid  of  any 
optical  devices.  Although  theoretical  studies  have  been  made  of  this 
phenomenon  and  of  the  intermediate  images  that  are  formed  ("Fresnel 
images")  it  has  not  been  fully  exploited  for  the  design  of  high-quality, 
low-cost  optical  instruments.  This  report  is  concerned  with  the  proper 
modifications  of  this  effect.  Two  new  shearing  interferometers,  an 
autocollimator,  a  Fourier  analyzer  and  a  spectrometric  imager  using  this 
effect  are  developed  in  this  report. 

Following  an  introduction  and  a  brief  review  of  the  pertinent 
details  of  the  Talbot  phenomenon  theory,  the  Talbot  shearing  interfero¬ 
meters  are  discussed.  The  theory  and  experimental  evidence  are  presented. 
Multiple-shearing  interferences  are  obtained  that  can  be  reduced  to 
triple-shearing  or  double-shearing  interferences  by  the  addition  of  simple 
spatial  filtering.  When  the  shear  is  less  than  the  width  of  the  details 
in  the  object,  these  interferences  become  either  the  second  or  first 
derivative  of  the  object  under  test,  respectively.  Either  lateral  or 
constant  radial  shear  can  be  introduced  by  choosing  Ronchi  rulings  or 
circular  gratings.  Thus  both  lateral  and  radial  derivatives  are  easily 
obtained.  If  white  light  is  used  as  a  source,  color  fringes  or  high 
contrast  are  observed. 

The  spectrometric  imager  follows  and  is  an  instrument  that  encodes 
the  brightness  distribution  of  a  scene  S  as  a  function  of  wavelength  and 
angular  position  coordinates  (a, 3)  by  a  single  detector  which  responds  to 
the  light  intensity  as  two  gratings  are  displaced.  The  situation  for 

vi 


encoding  only  the  wavelength  information  S(X)  conceived  by  Lohmann  with 
spectra  by  Klages  is  briefly  reviewed.  A  Fourier  analyzer  is  discussed. 
The  function  to  be  analyzed  is  cut  out  of  an  opaque  material  and  is 
incoherently  and  uniformly  illuminated.  The  detected  signal  is  its  cosine 
transformation  as  a  function  of  grating  separation.  This  is  followed  by 
the  theory  for  an  instrument  that  encodes  only  S(a,6).  The  encoding  for 
a  restricted  region  of  (a, 3)  is  a  cosine  transformation,  and  therefore 
another  cosine  transformation  recovers  the  information.  Since  only  one 
point  in  the  two  dimensional  Fourier  cosine  transformation  space  can  be 
encoded  at  any  given  time,  the  Fourier  space  must  be  sampled.  Three 
different  scanning  systems  are  explored.  Limitations  of  our  instrument 
in  encoding  S(a,3>X)  are  discussed.  The  additional  modification  required 
to  encode  this  most  general  situation  and  the  attendant  loss  of  light 
energy  at  the  detector  are  also  discussed. 

A  separate  chapter  concludes  the  body  of  the  report  by  summarizing 
the  work  and  offering  remarks  for  future  projects.  Three  papers  are 
included  in  the  report.  Two  of  them  are  our  early  results  on  the  shearing 
interferometers  and  contain  interf erograms  of  good  quality.  The  third 
paper  is  about  an  autocollimator  which  is  a  special  case  of  the  Talbot 
interferometer.  Here  both  grating  rulings  are  aligned  in  the  same 
direction,  and  moir^  fringes  are  observed  in  the  image  plane  whenever 
the  two  grating  periods  are  unequal.  These  fringes  indicate  that  a 
spherical  wave  is  at  the  first  grating.  By  adjusting  the  collimating 
lens,  this  spherical  wave  can  be  made  plane  through  observing  the  moir§ 
fringes  at  the  second  grating.  A  uniformly  bright  (or  dark)  field  indi¬ 
cates  that  the  system  is  aligned  within  the  accuracy  limitations  of  the 
interferometer . 
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CHAPTER  I 


INTRODUCTION 

A  rather  remarkable  optical  phenomenon  was  reported  by  H.  F, 
Talbot^  over  one  hundred  and  thirty  years  ago.  It  is  the  goal  of  this 
report  to  use  this  effect  for  several  new  metrological  purposes. 

Briefly  stated,  this  phenomenon  that  we  call  the  "Talbot  effect" 
concerns  the  formation  of  images  when  a  certain  class  of  structures  like 
a  grating  are  illuminated  by  spatially  coherent  light.  In  other  words, 
the  Talbot  effect  is  a  method  for  forming  images  of  these  structures 
without  any  lenses  or  mirrors,  simply  by  illuminating  them  with  spa¬ 
tially  coherent  light.  The  first  object  was  a  grating  of  equidistant 

lines,  and  subsequent  research  has  enlarged  the  class  of  objects  to 
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include  a  number  of  periodic  structures  .  Montgomery  determined  the 
total  class  of  objects,  and  showed  that  the  two  dimensional  Fourier 
transform  of  these  self-imaging  objects  could  exist  only  on  a  set  of 
concentric  circles.  The  self-images  repeat  themselves  at  intervals 
deduced  by  Lord  Rayleigh^  as  2d  /A  (where  d  is  the  grating  period,  X 
is  the  wavelength  of  the  light)  for  plane  wave  illumination.  Therefore, 
longer  wavelength  illumination  produces  closer  spaced  self-images  than 
shorter  wavelengths.  Therefore,  when  white  light  is  used  as  a  source, 
color  images  of  high  contrast  are  observed. 

The  wavelength  dependency  of  the  intervals  between  self-images 
was  used  to  develop  a  spectrometer  by  A.  Lohmann^.  The  recorded  inter- 
ferogram  is  a  cosine  transformation  of  the  spectrum  and  is  the  first 
grating  spectrometer  that  records  all  wavelengths  in  the  spectrum 
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simultaneously  •  The  Talbot  effect  has  been  employed  in  connection 

with  masks  for  integrated  circuits  that  contain  hundreds  of  identical 

shapes  in  a  periodic  array  .  When  a  section  of  the  mask  is  damaged 

through  mishandling  or  wear,  the  original  pattern  can  be  reconstructed 

by  illuminating  the  damaged  mask  with  a  plane  monochromatic  wave.  In  a 

distant  self-image  plane  the  restored  image  appears.  Thus  far  only 

simple  geometries  have  been  successfully  tested.  The  high  contrast 

color  self-images  can  be  used  artistically.  By  placing  a  second  grating 

in  the  field  and  slightly  tilting  it,  moire  fringes  are  formed  between 

the  self-images  and  the  grating.  The  location  of  the  grating  selects 

the  color  combination.  A  third  grating  adds  more  fringe  patterns  and 

additional  colors.  A  similar  setup  was  used  by  E.  Lau  in  the  design  of 

8 

table  cloths  . 

The  useful  aspects  of  the  Talbot  effect  in  the  design  of  new 

instruments  follow.  Poor  quality  gratings  can  be  used  because  the  self- 

images  emphasize  the  periodic  structure  of  the  grating,  and  therefore 

local  defects  in  it  are  minimized.  Furthermore  the  location  of  the  self- 

images  are  at  intervals  that  are  a  function  of  the  grating  period,  and 

misalignment  of  the  grating  reduces  the  effective  grating  period  by  a 

cosine  factor.  In  addition  the  self-images  occur  over  distances  that 

are  easily  managed  by  the  experimenter.  The  self-image  interval  for 

-3 

d  =  0.1mm  and  X  =  2/3  x  10  mm  is  30  mm.  Consequently  positioning  of 
components  can  be  done  without  the  aid  of  micropositioners. 

Although  some  work  has  been  done  using  the  Talbot  effect,  it  has 
not  been  fully  exploited  in  the  design  of  rugged,  high  quality,  low  cost 
instruments.  This  report  is  concerned  with  the  proper  modification 
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of  this  effect  in  the  design  of  such  instruments.  Two  new  shearing 
interferometers,  an  autocollimator,  a  Fourier  analyzer,  an  image  scanner 
and  a  spectrometric  imager  are  developed.  These  instruments  are  dis¬ 
cussed  within  these  pages.  Each  chapter  has  appropriate  introductory 
remarks.  Chapter  II  is  a  review  of  the  Talbot  effect.  Only  relevant 
aspects  of  the  theory  are  covered.  Both  the  lateral  and  constant  radial 
shearing  interferometers  are  discussed  in  Chapter  III.  Ronchi  rulings 
give  a  lateral  shear,  and  circular  gratings  produce  a  constant  radial 
shear.  Therefore  either  one  dimensional  or  axially  s}anmetric  optical 
systems  can  be  tested.  With  simple  spatial  filters  either  double  or 
triple  shearing  Interferences  are  obtained  that  become  the  first  or 
second  derivative  respectively  when  the  shear  is  less  than  the  width  of 
the  tiniest  details  in  the  test  object.  The  theory  with  experimental 
evidence  is  presented.  In  Paper  1  and  Paper  2  are  reprints  of  our  pre¬ 
liminary  work  with  these  Interferometers  and  contain  interf erograms  of 
high  quality.  Chapter  IV  contains  theoretical  studies  of  the  Fourier 
analyzer,  the  image  scanner  and  the  spectrometric  imager.  The  brightness 
distribution  S  as  a  function  of  its  angular  coordinates  01,3  and  its 
wavelength  X  is  encoded  onto  a  single  detector  by  the  spectrometric 
imager.  The  image  scanner  encodes  only  S(ci,3),  and  the  Fourier  analyzer 
is  a  special  case  of  the  image  scanner.  Chapter  V  concerns  our  con¬ 
clusions  about  these  Instruments  and  remarks  about  future  work.  The 
autocollimator  is  covered  in  Paper  3. 

These  instruments  based  on  the  Talbot  effect  have  a  lot  in  common 
and  yet  differ  in  ways  that  allow  special  measurements  to  be  made.  All  of 
the  instruments  use  two  gratings, G1  and  G2,  where  self-images  of 
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the  first  grating  are  detected  by  a  second.  But  they  differ  in  the  way 
that  the  optical  light  signal  is  processed.  Both  the  shearing  inter¬ 
ferometers  and  autocollimator  use  an  area  detector  like  the  eye  or  photo¬ 
graphic  film,  and  their  gratings  are  kept  at  a  fixed  distance,  whereas 
the  remaining  instruments  use  a  combination  of  grating  motions  to  mod¬ 
ulate  the  light  intensity  onto  a  single  detector.  The  motions  used  are 
moving  G2  laterally  and  longitudinally  and  rotating  both  gratings  in 
unison.  The  latter  motion  is  required  when  two  dimensional  spatial 
information  is  desired.  The  lateral  motion  makes  it  possible  to  select 
either  one  of  two  terms  in  the  point  spread  response  that  contains  the 
wavelength  or  spatial  information.  A  compilation  of  the  physical  aspects 
of  our  Talbot  instruments  follow  with  a  summary  in  Table  I.  The  spec¬ 
trometer  is  added  to  complete  the  list,  and  diagrams  of  the  instruments 
can  be  found  in  the  appropriate  chapters. 

Shearing  interferometer:  Collimated  light  impinges  upon  G1  to 
form  the  self-images.  A  second  grating  is  situated  in  one  of  the  self- 
image  planes.  Placing  a  test  object  between  the  two  gratings  distorts 
the  self-images  at  G2.  The  moire  fringes  formed  by  the  distorted  self- 
images  and  G2  are  observed  by  an  area  detector. 

Autocollimator:  This  instrument  is  identical  to  the  above;  the 
test  object  is  the  curvature  of  the  wavefront  at  G1  caused  by  defocusing 
the  collimator  lens.  The  moire  fringes  measure  the  wavefront  curvature 
which  is  related  to  the  amount  of  defocusing. 

Spectrometer:  The  source  to  be  measured  is  collimated  and  self- 

images  of  G1  are  formed.  G2  moves  laterally  and  the  instrument’s 

* 

response  is  filtered  for  the  (a±X/d)  dependent  term  with  a  =  0.  The 
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light  is  collected  onto  a  single  detector  which  records  the  light  intens¬ 
ity  as  G2  is  moved  longitudinally. 

Fourier  Analyzer:  A  function  mask  is  cut  out  of  opaque  material, 
incoherently  illuminated  and  placed  in  the  front  focal  plane  of  the 
collimating  lens.  G2  is  moved  laterally  so  that  the a-dependent  term  of 
the  response  is  selected.  The  light  is  collected  onto  a  single  detector, 
and  the  intensity  is  recorded  as  G2  moves  longitudinally. 

Image  scanner:  The  source  to  be  measured  is  collimated  and  waves 
emitted  from  the  collimator  impinge  upon  Gl.  G2  is  moved  laterally  so 
that  the  a-dependent  term  can  be  selected  from  the  point  spread  response. 
The  light  is  collected,  and  the  intensity  is  recorded  as  a  function  of 
the  grating  separation.  In  addition  both  gratings  are  rotated  so  that 
tfte  two  dimensional  scene  can  be  scanned. 

Spectrometric  imager:  This  instrument  has  the  same  setup  as  the 
image  scanner  except  that  the  (a  +  X/d)  term  is  selected. 

The  foregoing  has  been  summarized  in  Table  I.  The  terms  x(t), 
z(t)  and  <f>  (t)  refer  to  the  grating  motions:  G2  lateral,  G2  longitudinal 
and  both  gratings  rotational.  The  a  and  (a  ±  X/d)  columns  stand  for  the 
selected  term  of  the  point  spread  response.  The  type  of  detector  used 
is  indicated  by  the  area  and  single  columns. 
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TABLE  I.  A  Comparison  of  the  Talbot  Instruments 


Instruments 

area 

single 

x(t) 

z(t) 

a 

a±X/d 

Autocollimator 

XX 

— 

— 

— 

— 

— 

— 

Shearing  Interferometer 

XX 

— 

— 

— 

— 

— 

— 

Spectrometer 

— 

XX 

XX 

XX 

— 

— 

XX 

Fourier  Analyzer 

— 

XX 

XX 

XX 

— 

XX 

— 

Image  Scanner 

— 

XX 

XX 

XX 

XX 

XX 

— 

Spectrometric  Imager 

— 

XX 

XX 

XX 

XX 

«... 

XX 

XX  marks  the  features  found  in  the  instrument. 
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CHAPTER  II 
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REVIEW  OF  THE  TALBOT  EFFECT 
INTRODUCTION 

In  1836  H.  F.  Talbot,  observing  with  a  powerful  lens  the  region 
behind  a  grating  placed  ten  to  twenty  feet  from  a  radiant  point  of  in¬ 
tense  solar  light,  discovered  a  series  of  images  of  the  grating.  In 
other  words,  the  diffracted  light  formed  images  of  the  grating  without 
the  use  of  any  optical  components.  We  call  these  images  ”Talbot  images," 
which  are  also  known  in  the  literature  as  "self  images"  and  "Fourier 
images . 

Lord  Rayleigh,"^  in  1881,  deduced  that  for  a  grating  of  period  d, 

illuminated  by  a  monochromatic  plane  wave  of  wavelength  A,  a  series  of 

Talbot  images  are  formed  at  intervals  given  by  =  A/ (1-  A-A^/d^  )  ^ 

2d^/A  when  A/d  <<  1.  Recently,  a  theoretical  treatment  with  experimental 

evidence  of  Talbot  images  for  structures  other  than  a  linear  grating  was 

2 

reported  by  Cowley  and  Moodie.  The  existence  of  many  other  images 

which  occur  in  the  region  of  Fresnel  diffraction  was  also  discussed. 

These  "Fresnel  images"  have  been  theoretically  explained  by  Winthrop  and 
3 

Worthington.  The  complete  set  of  objects  that  self-image  in  a  plane 

L 

when  illuminated  by  a  plane  monochromatic  wave  was  derived  by  Montgomery. 
Only  those  aspects  of  the  theory  that  are  required  to  understand  this 
dissertation  will  be  presented  in  the  following  section. 


Theory 

Let  a  monochromatic  plane  wave  be  incident  upon  a  grating  as  shown 
in  figure  1.  We  chose  the  mathematically  simple  grating  transmission 
function  g(x)  =  1  +  cos(27Tx/d),  so  that  the  relevant  parts  of  the  theory 
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G1 


Figure  1.  A  grating  illuminated  with  a  plane  wave.  The  observer  sees 
the  ”Talbot  image”. 

will  be  enhanced.  The  complex  wavefield  just  behind  the  grating  is  given 
by  the  product  of  the  grating  transmission  function  and  the  wavefield 
just  prior  to  the  grating.  The  complex  wavefield  behind  the  grating  can 
be  described  by  the  Raleigh-Sommerf ield-Debye  formulation  for  diffraction 
of  plane  waves.  We  have 


u(x,y,z) 

where 

u(v,y;0''’) 


+  vx+yy]  I  dvdy 


u(x,y,0  )  exp|-2Tri(vx+yy)|  dxdy. 


and  the  integrations  are  over  the  interval  +  «)  in  this  chapter. 

Perfixrming  the  Fourier  integration  yields 

u(v,y;0‘'’)  =  6(v,y)  +  1/2  6(v-l/2,y)  +  1/2  6(v+l/2,y), 
and  placing  this  into  our  diffraction  formula  yields 

u(x,y,2)  =  exp{ik.2}  +  cos(2TTx/d)  expjikz  A-X^/d^  }  (1) 
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This  equation  can  be  interpreted  as  the  sum  of  two  interfering 
waves  traveling  in  the  +z  direction.  The  first  wave  is  an  ordinary 
plane  wave  with  wavelength  X  and  phase  velocity  c,  and  the  latter  wave 
has  a  wavelength  =  X//l-X^/d^  and  phase  velocity  v^  =  c//l-X^/d^.  We 
can  understand  the  Talbot  effect  as  the  interference  between  an  ordinary 
plane  wave  and  a  modulated  plane  wave.  Since  their  phase  velocities  are 
different,  they  must  travel  a  distance  before  they  are  again  in  phase. 

We  can  rewrite  equation  (1)  in  a  way  that  will  best  illustrate 
this  explanation: 

u(x,y,z)  =  exp{ikz}  [l  +  cos(2irx/d)  exp|ikz(l  -  /1-X^/ d^)| ] 

The  overall  phase  factor  exp  {ikz}  can  be  eliminated  from  further  con¬ 
sideration  since  the  final  observed  quantity  is  intensity.  Images  of 
the  grating  will  occur  whenever  the  exponential  factor  of  the  modulated 
wave  is  1 ;  and  this  occurs  at  multiples  of  the  distance 

z  =  x/(l  -  2d^/X  for  X/d  «  1. 

A  typical  example  for  z^  is:  d  =  0.1  mm,  X  =  2/3  x  10  mm;  z^  —  30  mm. 

Since  the  distance  for  Talbot  imaging  is  X  dependent,  this  effect  can  be 

6,7 

used  to  make  a  spectrometer. 

In  practice  gratings  are  of  limited  size,  and  the  diffracted 
orders  "walk-off"  the  zeroth  order  as  shown  in  figure  2.  The  Talbot 
effect  occurs  within  the  region  where  the  diffracted  orders  interfere. 
These  waves  overlap  in  the  shaded  region  of  figure  2,  which  forms  a 
wedge  whose  width  decreases  to  zero  at  z^^  =  Bd/2X,  where  B  is  the 
grating  width,  and  is  the  maximum  distance  from  the  grating  over 

which  the  Talbot  effect  occurs. 
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Figure  2.  The  first  diffracted  orders  "walkoff"  the  zeroth  order  with 
increasing  z.  The  maximum  distance  for  the  Talbot  effect 
is  marked  z^. 
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CHAPTER  III 


THE  TALBOT  SHEARING  INTERFEROMETERS 

INTRODUCTION 

Shearing  interferometry  has  been  used  to  test  optical  instruments 
and  components,  to  study  physico-chemical  phenomena  in  liquids,  and  to 
evaluate  optical  transfer  functions,  i.e,  the  contrast  of  an  image 
plotted  against  its  spatial  frequency.^  In  this  type  of  interferometry 
the  wavefront  under  test  is  duplicated  and  displaced.  They  are  then 
made  to  interfere  with  each  other,  yielding  a  comparison  of  the  wave- 
front  with  a  sheared  image  of  itself.  These  interferometers  typically 
use  a  single  beam-dividing  element  that  affords  a  simplicity  over  other 
interferometers.  There  are  two  ways  of  observing  the  sheared  wavefronts 
The  first  method  is  to  have  the  shear  sufficiently  large  that  each  of 
the  wavefronts  interferes  with  the  uniform  background  of  the  other.  This 
is  called  "total  shear"  and  the  fringes  observed  are  lines  of  equal 
optical  paths.  In  the  second  method  the  shear  is  kept  smaller  than  the 
width  of  the  details  that  we  wish  to  observe,  and  the  fringes  are  pro¬ 
portional  to  the  gradient  of  the  wavefront  in  the  direction  of  the  shear 

This  method  is  known  as  "differential  shear."  The  earliest  shearing 

2 

interferometer  introduced  a  lateral  displacement  between  the  waves,  and 

the  most  recent  class  of  shearing  interferometers  cause  a  constant 

3 

radial  shift  between  the  wavefronts.  In  the  former  case  we  obtain  the 
lateral  gradient  of  the  wavefront,  and  in  the  latter  the  radial  gradient 
is  observed.  It  is  desirable  to  have  both  types  of  interferometers 
since  geometries  that  require  testing  are  often  either  one-dimensional 
or  circular. 
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Recently,  a  new  shearing  interferometer  based  on  the  Talbot  effect 

A 

was  introduced.  The  basic  setup  is  shown  in  fig.  4.  An  image  of  the 
grating  is  formed  at  multiples  of  a  distance  =  2d^/A  (where  d  is  the 
grating  period)  behind  the  grating  that  has  been  illuminated  by  a  plane 
monochromatic  wave.  The  space  between  the  two  gratings  is  empty.  In 
other  words,  the  Talbot  effect  is  a  method  of  forming  images  of  a  peri¬ 
odic  object  without  any  lenses  or  mirrors.  When  properly  modified,  the 
Talbot  effect  affords  a  unique  method  for  quantitatively  observing  phase 
objects  through  the  use  of  two  gratings,  one  of  which  is  self-imaged 
onto  the  second.  Since  the  grating  shears  the  wavefront  in  a  direction 
perpendicular  to  its  rulings,  a  Ronchi  ruling  laterally  shears  the  wave- 
front  and  a  circular  grating  gives  a  constant  radial  shear  to  the  wave- 
front.  Therefore,  by  a  suitable  choice  of  gratings,  optical  instruments 
with  either  one-dimensional  or  circular  geometries  can  be  tested  simply. 

In  this  paper  the  complete  theory  and  experimental  evidence  for  both 
types  of  shearing  interferometers  will  be  presented,  A  complete  theory 
requires  scalar  diffraction  theory  to  explain  the  self-imaging  effect, 
the  spatial  filtering  capability,  and  the  high-contrast  color  fringes 
that  are  observed  when  white  light  is  used  as  a  source.  However,  a  ray 
explanation  like  that  of  Nishijima  and  Oster^  will  be  presented  first 
to  acquaint  the  reader  with  the  operation  of  this  instrument.  The  ray 
explanation  is  valid  only  for  distances  from  the  first  grating  that  are 
about  one-eight  of  the  Talbot  distance  z^,  which  greatly  limits  the  sensi¬ 
tivity  by  orders  of  magnitude.^  Thus  the  ray  explanation  will  be  fol¬ 
lowed  by  the  wave  theory  in  order  to  understand  the  full  capabilities 
of  this  instrument. 
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A  Qualitative  Ray  Explanation 


A  Ronchi  ruling  is  illuminated  with  collimated  light  as  shown  in 
Figure  3.  The  diffracted  waves  are  modulated  and  form  "Fresnel  Images." 
At  integer  multiples  of  the  distance  z^«  2d^/X,  where  d  is  the  grating 
period,  an  image  of  the  Ronchi  ruling  is  formed.  If  a  second  grating  is 
placed  at  one  of  these  image  planes  in  anti-position  to  the  image,  no 
light  will  be  transmitted  beyond  the  second  grating.  Hence  a  ray  that 
passes  through  a  slit  of  the  first  grating  is  blocked  by  the  second 
grating  as  shown  in  Figure  3.  A  prism  of  wedge  angle  a.  is  introduced  at 
a  distance  z  from  the  second  grating  and  normal  to  the  Incident  rays. 

The  rays  passing  through  the  prism  will  be  bent  through  an  angle  e.  Thus 
the  intensity  of  the  light  as  seen  by  the  observer  varies  with  the  angle 

e,  with  peak  intensities  given  by  **  tan  Ejh  “  (m  +  1/2)  d/z,  and  dark 

fringes  by  tan  =  md/z.  The  angle  e  is  related  to  the  wedge  angle 

a  under  small  angle  approximation  by  e  =  (n  —  l)a.  The  wedge  angle  ct  is 
for  small  angles  determined  by  the  ratio  a  ^  tan  a  =  At/Ax,  where  At  - 
prism  thickness.  Ax  =  distance  from  apex.  A  phase  object  can  be  thought 
to  be  made  up  of  prisms,  a  fact  which  is  often  used  to  explain  the  action 
of  a  lens,  so  that  for  any  point  (x,y)  on  the  object,  the  prism  angle  a 

is  given  by  ct(x,y)  ^  8t(x,y)/9x,  which  in  turn  is  related  to  e  by  e(x,y) 

=  (n-1)  3t(x,y)/9x.  The  derivative  9t(x,y)/9y  is  of  no  concern  here  since 
the  grating  bars  are  parallel  to  the  y-axis. 

•  Thus  a  dark  fringe  structure  is  obtained  whenever  e(x,y)  meets  the 
dark  fringe  conditions.  These  fringes  indicate  lines  of  equal  9t/9x, 
which  is  the  x-component  of  the  surface  gradient  at  point  (x,y) .  Where 
the  lines  are  close  together,  the  slope  9t/9x  changes  rapidly.  Hence 
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the  fringe  density  or  fringe  frequency  Is  proportional  to  the  x-component 
of  the  curvature  8^t/3x^.  The  maximum  detectable  curvature  Is  determined 
by  the  minimum  detectable  fringe  separation.  In  a  sense  what  Is  observed 
Is  the  second  derivative  of  the  phase  object. 

There  are  two  problems  of  fringe  detection.  The  first  occurs  when 
the  fringes  become  too  fine  for  the  eye  to  resolve  them.  This  Is  only 
a  practical  limitation  and  a  magnifying  glass  would  extend  the  sensiti¬ 
vity.  But  the  second  problem  Is  more  fundamental.  The  field  of  view  Is 
finite,  say  B  =  Nd,  where  N  Is  the  number  of  lines  In  the  grating.  For 
an  object  whose  curvature  Is  at  the  minimum  detectable  sensitivity,  Its 
curvature  must  be  large  enough  for  at  least  one  full  (or  maybe  one-half) 

fringe  across  the  field.  Hence  [(9t/3x)  -  (3t/3x)  .  1  (n-1)  >;  d/2z. 

max  min 

If  (3t/9x)  Is  at  X  =  +Nd/2  and  (9t/9x)  .  Is  at  x  =  -Nd/2,  then 
max  '  '  min  ’ 

3t(x+Nd/2,y)/3x  -  9t(x-Nd/2,y)3x  ««  9^t/3x^  Nd. 

d  /  2  X 

Hence  |d^t/9x^|  ^  z(n"l)Nd  ”  yNz(n’-l)  ’  maximum  Talbot  distance  is 

given  by  z  =  Nd^/ZX.  Using  this  as  z  In  the  above  equation  yields 
max 


Is^t/ 


(Nd) ^ (n-1)  B^(n^l)  * 

For  example  If  X  =  5  x  10“^*  mm,  n  -  1  =  1/2,  B  =  10^ 


|92t/9x2|  -  10“^ 


mm 


■  1  _ 


=  lAlO  km).  We  note  that  z 


mm,  we  obtain 
can  be  Increased  with 


max 

the  addition  of  plane  mirrors  that  are  positioned  tangent  to  the  zeroth 
order  beam  which  effectively  Increases  the  grating  size  due  to  reflection. 

As  an  example  of  the  system,  a  perfect  thin  lens  Is  considered. 

The  phase  of  this  lens  Is  given  by  (p(x,y)  =  (-tr/Xf)  (x^  +  y^),  where  f  Is 
the  focal  length  of  the  lens,  and  (x,y)  Its  coordinates.  The  relative 
phase  of  any  object  Is  given  by  (p(x,y)  =  (2t:/X)  (n-1)  t(x,y),  and  there- 
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fore  t(x,y)  =  -(x^  +  y^)  /  [2f  (n-1)  ]  .  Thus  9t(x,y)/8x  =  -x[(n-l)f],  which 
when  substituted  into  equation  (1)  gives  e(x,y)  -  -x/f.  This  is  a 
straight  line  that,  when  compared  with  conditions  for  dark  fringes,  will 
yield  equally  spaced  fringes  of  fd/z.  An  imperfect  lens  will  not  have 
straight  equally  spaced  fringes,  and  thus  its  quality  can  be  determined. 
It  is  worth  noting  that  this  system  likewise  tests  the  collimator  objec¬ 
tive  and  can  serve  as  a  good  method  for  testing  lenses. 

The  Lateral  Shearing  Interferometer 


Discussion 

In  this  section  a  rigorous  derivation  based  on  scalar  diffraction 
theory  is  presented.  The  result  will  be  an  interpretation  of  the  image 
as  two,  three, or  many  shift ed-obj ect  wavefronts,  sometimes  tilted.  Ac¬ 
cordingly  we  will  use  the  terms  **shearing  interferometry,”  as  is  common 
when  two  shif ted-object  wavefronts  interact.  In  generalization  thereof 
we  will  introduce  the  terms  ^triple  shearing  interferometry”  and  ”mul- 
tiple  shearing  interferometry.”  In  some  instances  these  shearing  inter¬ 
ferences  will  represent  the  first  or  the  second  derivative  of  the  object. 
The  instrument  is  diagrammed  in  figure  4.  The  spatial  coherence  require¬ 
ment  is  derived  in  Appendix  1. 

Theory 

Consider  a  monochromatic  plane  wave  incident  on  grating  G1  in  the 
plane  at  -z^.  The  field  behind  this  grating  is 

n=-oo 

where  the  periodic  grating  is  considered  infinite  in  extent  and  of  period 
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d,  and  is  expressed  as  a  Fourier  exponential  series j  k  —  2'ir/X  where  X 

I  8  , 

wavelength  of  light.  Following  a  line  of  reasoning  like  Edgar  s,  the 
incident  plane  wave  is  by  action  of  the  grating  replaced  by  a  set  of 
plane  wavefronts  whose  x— direction  cosines  are  nX/d  and  whose  field 
strength  is  proportional  to  The  propagating  field  behind  G1  is 


.  ,  -ikzi\  {ik[nxX/d  +  (zj^+z)  /l-(nX/d)^] } 

u(x,y,z)  =  e  e 

(n) 


which  just  prior  to  the  object  plane  z  =  0-  reduces  to 


u(x,y,0-) 


{ik[nxX/d  +  z^ /l-(nX/d)^] } 
e  ^ 


Applying  the  Kirchhoff  boundary  conditions  to  the  field  as  it  passes 

through  the  object,  we  have 

-Ikzi  \  '  \  {ik[nxX/d  +  z, /l-(nX/d)  ^]  } 

u(x,y,0+)  =  e  12^C^uo(x,y)  e  1 

(n) 


UQ(x,y)  is  tli6  two— diin6nsiona.l  object  transtnittance  function. 

Each  of  the  plane  waves  is  diffracted  as  it  passes  through  the  ob¬ 
ject.  At  this  point,  the  concept  of  the  angular  spectrum  is  very  useful. 
In  this  formulation  the  diffraction  phenomenon  is  a  multiplicative  quad¬ 
ratic  phase  factor  exp  [ikz/l-X^  (v^+  y^)  ]  the  Fourier  domain  which  in¬ 
creases  with  the  propagation  distance  z.  Let  u(v,y;  0“h)  be  the  angular 
spectrum  of  equation  (2)  defined  by 


u(v,v;  0+)  =JJu(x,y,  0+)  e  ^iriCvx+yy) 
where  the  integration  in  this  section  will  be  over  the  interval 
when  unspecified,  and  (v,y)  are  the  Fourier  spatial  frequency  components. 
The  angular  spectrum  for  the  plane  just  preceding  G2  (z  =  Z2-)  is 
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u(v,ii;  z  -)  =  u(v,y;  0+)  exp 


{  ikz2'^l-X^(v^+y^) } 


^^-ikzi^  C^uo(v-n/d,y)  ^WziA-inX/dY^  + 
(n) 

Z2A~X^(v^+y^)  ] } 

where  uo(v,ij)  is  the  angular  spectrum  of  the  object.  The  Kirchhoff 
boundary  conditions  are  applied  to  the  second  grating  at  z  =  Z2 


u(x,y,Z2+)  =  u(x,y,Z2 


-)E 

(m) 


C  e 
m 


iTTm  27riinx/d 


where  the  factor  accounts  for  the  fact  that  the  second  grating  is 

shifted  by  half  a  period  with  respect  to  Gl.  In  the  Fourier  domain  this 
becomes  a  convolution  which  results  in 


u(v,y;  Z2+)  =  e  ^  ^  ^  uq  (v-(n-Kn) /d,y) 

(n)  (m) 

ik{ z  1  /l  - (nX/ d)  ^  +  Z2  A-X^ [  (v-m/d)^+y^]  }}  (3) 

The  image  plane  is  conjugate  to  the  object  plane  z  =  0.  There¬ 
fore,  we  compute  the  field  virtually  back  to  the  z  =  0  plane  in  order 
to  find  the  resultant  field  in  the  plane  of  observation,  and  give  the 
field  a  new  symbol  v(x,y)  with  v(v,y)  as  its  Fourier  spectrum 
v(v,y)  =  u(v,y;  Z2+)  e{-ikz2 A-X^(v^+y^)  }  ^ 

Substituting  equation  (3)  into  this  and  simplifying  by  using  the  first 
two  terms  in  the  Taylor  series  expansion  of  the  square  roots  in  the  ex- 
ponential  terms  yields 


v(v,y) 


^  ^  ,  ^n^m®^™  uo(v-(n+m)/d,y) 

(n)  (m) 

^{ik[z2A2vm/d  -  Zj(nx/d)2  +  Z2(mx/d)2]}. 
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But  as  explained  in  the  first  section,  the  two  gratings  are  separated 
by  some  multiple  M  of  the  Talbot  distance  z^.  Thus  putting  Z2  +  Zj^  = 
M  z^  «  2Md^/X  into  the  above  equation  reduces  it  to 


v(v 


,y)  =  /  J  /  i  uo(v-(n+m)/d,y) 


(n)  (m) 


g{2TTi[Az2vm/d-Z2X(m2-n2)  /2d^]}.  (4) 

Taking  the  Fourier  transform  of  this  gives  the  resultant  field  in  the 


image  plane. 


v(x,y) 


(n)  (m) 


C  C  Ur,  (x+mz2X/d,y) 
n  m  ^ 

{2iri[m/2  +  (m+n)2z2/z  +  (n+m)x/d]}. 

e 

(5) 


Tilted  Ronchi  Rulings 

Thus  far  our  analysis  indicates  that  the  system  cannot  distin¬ 
guish  a  positive  from  a  negative  lens,  since  it  is  insensitive  to  the 
sign  of  the  phase  derivative.  By  rotating  the  grating  of  Gl,  its  self- 
images  will  correspondingly  rotate,  and  moire  fringes  between  G2  and 
these  self-images  will  be  observed. The  frequency  of  the  moire  fringes 
and  their  orientation  relative  to  the  gratings  are  dependent  upon  the 
angle  a  between  the  grating's  lines  and  upon  each  grating's  frequency.  Un¬ 
equal  grating  frequencies  cause  rotation  of  the  moire  fringes,  a  fact 
that  was  exploited  in  making  a  sensitive  autocollimator  (see  Papers  3) . 

In  our  autocollimator,  a  diverging  wave  at  Gl  causes  the  self-image  to 
be  magnified,  which  corresponds  to  a  decrease  in  its  frequency,  whereas 
a  converging  wave  increases  the  frequency.  In  both  cases  the  Moire 
fringes  rotate,  but  in  opposite  directions.  Thus,  by  the  rotation  of 
the  fringes,  the  sign  of  the  curvature  of  the  wave  at  Gl  could  be  deter- 
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Figure  5(a) . 


The  reciprocal  space  for  calculating  the  moir^  effect 
when  I I  >  |  V2  |  . 


Figure  5(b) .  The  reciprocal  space  for  calculating  the  moire  effect 
when  I vj I  <  |  V2  |  . 
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mined.  This  phenomenon  can  be  simply  understood  using  a  technique  by 

T,  li 

Rogers  . 

Rogers  describes  a  reciprocal  space  that  is  like  a  Fourier  space 
where  vectors  in  the  space  describe  the  frequencies  of  a  grating.  Each 
vector  is  weighted  by  its  associated  Fourier  coefficient.  The  length  of 
the  vector  corresponds  to  the  value  of  the  grating  frequency  and  its 
orientation  is  perpendicular  to  the  grating  lines.  Our  diagram  contains 
vectors  for  the  self-image  of  G1  and  G2.  We  consider  only  the  funda¬ 
mental  frequency  of  each  grating,  since  they  are  the  dominant  ones  in 
the  Moire  effect.  The  vector  difference  yields  the  Moire  frequency 
vector. 

We  rotate  G1  by  -q)/2  and  G2  by  -Kp/2.  Let  v^^  be  the  vector  for 

the  self-images  and  V2  for  G2.  For  a  plane  wave  vi  =  V2  and  we  have  Moire 

fringes  parallel  to  the  x-axis.  When  the  self-images  are  minified,  as 

when  G1  is  illuminated  with  converging  wave  vi  >  V2,  we  have  the  vector 

diagram  of  Fig.  5(a).  The  Moire  vector  v^  has  rotated  counter-closewise 

from  its  vertical  position  for  vi  =  V2.  On  the  other  hand,  with  magnified 

self-images,  v^  <  V2  and  the  vector  diagram  of  Fig.  5(b)  results.  In 

this  case,  v,  has  rotated  clockwise.  By  choosing  to  rotate  G1  by  +cp/2 
M 

and  G2  by  -<p/2,  the  above  situations  would  be  reversed  with  clockwise 
rotation  of  the  Moire  fringes,  indicating  a  converging  and  counter-clock¬ 
wise  rotation  indicating  a  diverging  wavefront. 

When  test  objects  are  placed  in  the  field,  the  same  results  apply. 
12 

Yokozeki  and  Suzuki  have  an  expression  for  the  locus  of  these  fringes 

for  the  Talbot  interferometer 

^  cos  g-l  ^  ^  3g(x,,y)  +  _4_  ^ 

^  sin  a  sin  a  9x  sin  a 
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where  a  is  the  angle  between  ruling  lines  defined  as  positive  if  G2  is 
oriented  clockwise  to  Gl,  Z  is  the  object  to  G2  distance,  g(x,y)  is  the 
distortion  of  the  zeroth  order  diffraction  wave  at  G2  and  m  is  an  integer 
This  expression  shows  that  the  straight-line  Moire  fringes  are  altered 
in  proportion  to  the  first  derivative  of  the  distortion.  The  direction 
that  the  fringes  rotate  is  given  by  the  sign  of  our  derivative.  As  an 
example,  consider  an  ideal  negative  lens.  The  partial  derivative  is 
proportional  to  x,  which  would  increase  the  value  of  the  slope  for  the 
straight-line  function  resulting  in  a  counter-clockwise  rotation.  A 
positive  lens  would  decrease  this  slope,  giving  the  fringes  an  opposite 
rotation. 

Filtering 

Thus  multiple^shearing  interferences  are  observed  in  the  image 
plane.  In  ordinary  shearing  interference  only  two  terms  are  present, 
u(x  +  Ax,y)  -  u(x,y).  In  order  to  reduce  the  large  number  of  terms  of 
the  double  series,  we  will  perform  certain  spatial  filtering  operations. 
The  field  expressed  by  equation  (4)  is  found  in  the  spatial  plane  located 
one  focal  length  behind  the  first  lens  of  the  telecentric  system.  In 
order  to  filter  just  the  n+m=0orn+m=l  term,  which  we  will  show 
leads  to  special  interferences,  we  must  assume  that  the  shifted  spectrums 
uo (^-(n+m) /d,y )  do  not  overlap.  This  is  so  if  the  bandwidth  restriction 
of  the  object  1av|  <l/2d,  i.e.  uo(v,p)  =  0  whenever  |v|  ^l/2d.  It  can 
also  be  shown  that  allowing  both  first  orders  to  pass  leads  us  to  the 
auto-correlation  of  UQ(x,y), 
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a.  Zero  Order  Filtering 

The  multiple  shearing  interferences  can  be  simplified  by  spatial 
filtering.  Our  first  consideration  is  the  case  of  zero  order  filtering, 
n  +  m  =  0.  Thus  eq.  (5)  reduces  to 

^m^-m®  UQ(x+mA,y), 

(m) 

where  A  =  shear  introduced  by  the  grating  to  the  wavefront 

under  test.  If  we  restrict  ourselves  to  only  considering  either  ampli¬ 
tude  or  phase  modulated  gratings  that  are  symmetrical  about  the  origin, 

then  the  C  's  are  real  and  C  =  C  .  Thus  we  can  rewrite  the  above 
m  m 

equation  as 

v(x,y)  =  cg.  u”o(x,y)  +  2  ^  [uo(x+mA,y)  + 

uo (x-mA,y)] |. 

These  multiple-shearing  Interferences  can  be  simplified  to  triple-shear- 

ing  interferences  if  terms  are  all  but  for  the  m  =  1  term 

negligibly  small.  By  using  a  special  grating  g(x)  =  1  +  cos(2Trx/d), 

triple-shearing  interferences  would  result  since  =  0  for  all  m>l. 

We  now  consider  the  use  of  an  inexpensive  Ronchi  ruling  as  a  grating. 

The  Fourier  coefficients  are  C  =  (1/2)  sinc(m/2);  all  even  C  -  0, 

m  ^ 

where  sinc(x)  =  sin(Trx) /ttx.  With  this  grating  the  multiple-shearing 
interferences  appear  as  triple-shearing  interferences  since  (C^/Cq)  = 

0  for  m  >  1.  Note  that  the  highest  term  neglected,  n  =  3,  is  just  one 
ninth  of  the  smallest  term  that  we  kept,  m  =  1.  Therefore,  our  equa¬ 
tion  reduces  to 
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v(x,y)  C§{uo(x,y)  -  (Ci/Co)^  [uo(x+A,y)  +Uo(x-A,y)]} 

Here  we  have  triple-shearing  interferences  in  the  image  plane  of  the 
instrument  due  to  zero  order  filtering.  These  triple-shearing  inter¬ 
ferences  become  the  second  derivative  of  the  wavefront  when  the  shear 
is  made  very  small.  To  show  this  we  expand  UQ(x±A,y)  by  its  Taylor  ser¬ 
ies,  and  substitute  it  into  the  above  equation  to  yield 

00 

v(x,y)  ««Uo(x,y)  [C^-2C2]  -  2C§  (x,y)  A^/p!  (6) 

p=2 

where  uo^\x,y)  =  3^uo(x,y)/9x^.  In  order  to  eliminate  the  unshifted 
image  from  the  field  of  view  we  must  have  C§  =  2cf.  This  can  be  done  by 
choosing  a  grating  whose  mark-to-space  ratio  is  1.3.  If  one  uses  the 
Ronchi  ruling,  this  term  is  still  negligible  as  compared  to  the  second 
term  of  eq.  (6) .  This  summation  term  can  be  shown  to  be  approximately 
the  second  derivative  by  using  a  theorem  by  Bernstein  which  states 
that  the  derivatives  of  a  bandlimited  function  are  bounded  by  the  maxi¬ 
mum  value  of  the  function  in  the  interval.  The  actual  manipulations  are 
done  in  Appendix  2,  the  outcome  of  which  states  that  the  maximum  shear 
A  must  be  less  than  one-sixth  of  the  width  of  the  smallest  details  in  the 
object  in  the  direction  of  shear.  This  in  turn  determines  the  distance 
Z2  that  the  object  must  be  placed  before  the  second  grating.  For  spatial 
filtering  we  require  that  Av  -  1/d,  where  the  most  stringent  condition 
is  the  equality.  This  condition  is  necessary  to  keep  the  shifted  spec- 
trums  in  the  Fourier  domain  of  our  telecentric  system  from  overlapping. 
Since  the  shear  A  =  <  6x/6  =  l/6Av  =  d/6,  we  have  Z2  <  d^/6X  = 

z^/12.  For  a  Ronchi  ruling  of  10  lines/mm,  this  means  that  the  object 
must  be  placed  within  5  mm  of  the  second  grating.  With  such  a  placement 
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of  the  object  and  a  ruling  with  mark-to-space  ratio  of  unity,  eq.  (6) 
becomes 

v(x,y)  »«(l/20)  uo(x,y)  -  (1/10) (zX/d)^  (x,y)/3x2, 

and  with  a  grating  of  mark-to-space  ratio  of  1.3,  we  obtain 

v(x,y)  ~  92uo(x,y)/8x2.  Thus  the  result  of  zeroth  order  fil¬ 
tering  is  the  second  derivative  of  the  wavefront  under  test. 

b.  First  Order  Filtering 

The  mathematics  for  this  case  follows  simply  from  the  previous 
section.  Here  we  filter  either  the  n  +  m  =  1  or  n  +  m  =  -1  term.  The 
only  difference  is  in  the  sign  of  a  phase  factor,  which  is  unimportant 
when  detecting  the  signal.  For  the  n  +  m  =  ■H,  the  field  in  the  image 
plane  becomes 

v(x,y)  =  CgCie^  l[x/d+z/zj]  -  ug (x+zX/d,y) } . 

This  surprising  result  occurs  since  the  even  coefficients  are  zero. 

Those  terms  which  contribute  under  the  n  +  m  =  1  condition  are  = 

C  C  ,  ,  implying  that  for  any  m  0,  one  coefficient  will  always  be  even 
m  m-l 

and  hence  equal  to  zero.  This  is  very  much  like  ordinary  shearing  inter¬ 
ferometry,  and  we  will  again  show  that  under  a  certain  condition  this  is 
approximately  the  derivative  of  the  object.  We  ignore  the  phase  since 
it  is  the  intensity  which  is  observed  in  the  image  plane,  and  deal  with 
the  bracketed  terms.  The  last  term  is  expanded  in  a  Taylor  series  around 

zX/d;  combined  with  ug(x,y)  it  gives 

00 

{...}  =  UQP\x,y)  (zX/d)^(l/p!)  »«Uq^^  (x,y)  (zX/d)  . 

P-1 

Using  the  same  kind  of  arguments  as  in  the  previous  section,  this  approxi- 
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mation  is  valid  if  the  shift  zX/d  -  I/SttAv  ^  (1/15) 6x.  Under  the  stric¬ 
test  condition  |Av|  =  1/d,  this  means  that  the  object  to  G2  distance  is 
z  5  z^/IOtt.  Therefore,  the  field  in  the  image  plane  for  first  order 
filtering  under  these  two  conditions,  whichever  may  apply,  is  given  by 
v(x,y)  (zX/2d7T)  (9uo(x,y)/3x)  , 

which  is  essentially  the  first  derivative  of  the  object.  We  ignore  the 
multiplicative  phase  factors  since  we  detect  intensity  and  note  that  the 
first  order  filtering  gives  us  essentially  the  first  derivatives  of  the 
object  v(x,y)  Suq  (x,y)/9x. 

Experiments 

The  Talbot  shearing  interferometer  was  assembled  on  a  laboratory 
bench  using  two  Ronchi  rulings  of  10  lines  per  mm  pitch  as  gratings. 

The  second  grating  was  mounted  in  a  holder  that  could  be  rotated  or 
shifted  laterally,  thereby  allowing  the  gratings  to  be  aligned  as  desired. 
Both  a  helium-neon  laser  and  a  tungsten  lamp  were  used  as  light  sources. 
Color  fringes  of  high  contrast  were  obtained  with  the  tungsten  lamp  but 
will  not  be  presented  because  of  the  cost  of  color  printing.  The  output 
of  a  laser  was  collimated  with  a  6  micron  pinhole  and  250  mm  focal  length 
objective  producing  a  50  mm  beam.  The  second  grating  was  located  in  one 
of  the  self-image  planes  of  the  first.  Objects  to  be  tested  were  in¬ 
serted  between  the  gratings  at  a  distance  from  the  second  grating  that 
introduced  the  desired  shear.  The  method  of  observing  interference 
fringes  was  to  project  the  object  with  its  fringe  structure  onto  a  screen. 
Thus,  the  fringes  overlap  the  corresponding  regions  of  the  object  to 
which  they  apply. 
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Figure  6.  A  lens  tested  in  the  lateral  shearing  interferometer  showing 
fringes  of  equal  phase  derivative.  Horizontal  moire  fringes 
appear  in  the  background. 


27 


Figure  7(a).  Zeroth  order  filtering  of  the  wavefront  that  has  passed 
through  the  medium  heated  by  a  candle  flame.  The  three 
candle  wicks  verify  that  triple  shearing  interferences 
are  observed. 


Figure  7(b).  First  order  filtering  of  the  same  object  as  in  (a), 
shearing  interferences  are  observed. 


Double 
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The  case  of  a  lens  was  discussed  in  the  opening  section.  Straight 
fringes  were  predicted  for  a  lens  with  only  quadratic  phase  terms; 
this  is  shown  in  Figure  6.  One  feature  not  discussed  previously  is  the 
addition  of  the  moire  fringes  in  the  background.  This  is  caused  by  a 
slight  rotation  of  the  second  grating  with  respect  to  the  self-image  of 
the  first.  This  in  turn  causes  the  fringes  to  rotate. 

Heated  gases  can  be  studied  by  this  technique  since  the  gas  will 
have  a  refractive  gradient.  A  candle  flame  was  placed  before  Gl.  We 
used  both  a  uniform  dark  field  and  a  Moire  pattern  for  our  background 
while  observing  the  nimbus  of  the  candle.  These  pictures  have  been  pub¬ 
lished  in  our  earlier  publication  in  Paper  1  figs.  3(a)  and  3(b).  The 
addition  of  the  Moire  fringes  enhances  our  ability  to  observe  the  fringe 
location  and  gives  us  the  sign  of  the  derivative. 

As  a  demonstration  of  spatial  filtering,  this  same  candle  flame 
was  spatially  filtered.  The  shear  was  made  excessively  large  so  that 
triple-  and  double-shearing  interferences  could  be  recognized  by  counting 
the  replication  of  some  unique  section  of  the  object  such  as  the  wick  of 
the  candle.  In  the  experiment  the  candle  flame  was  separated  from  the 
second  grating  by  90  cm,  causing  a  shear  of  about  5.6  mm.  Zeroth-order 
filtering  is  demonstrated  in  figure  7(a).  Three  wicks  are  clearly  ob¬ 
servable  indicating  triple-shearing  interferences.  The  two  sheared  wicks 
are  of  lesser  contrast  than  the  unsheared  ones,  as  predicted  by  the  theory. 
In  figure  7(b)  we  have  first-order  filtering  which  produces  double-shearing 
interferences.  These  two  candle  wicks  are  of  equal  contrast  as  expected 
from  our  theory. 
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Constant  Radial  Shearing  Interferometer 


Discussion 

In  our  previous  instrument  using  Ronchi  rulings,  the  derivative 
of  a  wavefront  perpendicular  to  the  ruling  lines  was  displayed  by  a  set 
of  interference  fringes.  Imagine  that  we  possess  a  "flexible”  Ronchi 
ruling.  Bending  this  flexible  grating  into  a  complete  circle  gives  us  a 
circular  grating  where  the  derivative  is  still  perpendicular  to  the  lines, 
but  the  lines  are  now  concentric  circles.  Therefore,  we  expect  that 
such  a  grating  in  a  Talbot  interferometer  setup  would  give  us  the  radial 
derivative. 

In  the  succeeding  sections  we  will  learn  that  for  experimental 
reasons  the  first  four  periods  of  a  grating  are  ignored  (r  ^  4a)  and  that 
self-images  of  a  circular  grating  exist  outside  of  a  region  defined  by  a 
cone  along  the  optical  axis.  The  radius  of  this  cone  increases  as  we 
move  farther  from  the  grating  (r  =  Az/a).  Combining  both  restrictions 
gives  us  r  <  4a  +  z/a.  This  region  is  excluded  when  evaluating  the 
output  of  our  new  instrument. 

As  we  will  see,  the  constant  radial  shearing  interferometer,  like 
the  lateral  shearing  interferometer,  gives  color  fringes  of  high  contrast 
when  white  light  is  used  as  a  source,  and  yields  double  or  triple  shearing 
interferences  with  spatial  filtering.  These  filtered  interferences 
become  the  first  and  second  radial  derivative  of  the  wavefront,  respec¬ 
tively,  when  the  shear  is  less  than  the  width  of  the  object  details  of 
interest.  Experimental  evidence  will  be  presented. 
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Definition  of  a  Circular  Grating 


A  circular  grating  consists  of  a  number  of  concentric  circular 
rings  on  a  transparent  surf ace. Every  other  ring  is  opaque,  and  the 
ratio  of  opaque  rings  to  transparent  ones  is  chosen  by  the  experimenter 
in  a  fashion  that  will  be  discussed  under  spatial  filtering.  The  trans¬ 
mittance  of  the  grating  G(t)  shown  in  Figure  8  is  given  by 


G(r) 


2TTimr/a 

e  , 


4a  <  r  ^  A 


otherwise 


(7) 


where  the  grating  is  represented  by  its  Fourier  exponential  series  »  ^m 
ctsinc  ma,  r  is  the  radial  location  of  any  point  on  the  grating,  a  is  the 
periodicity  of  the  structure,  ota  is  the  dimension  of  the  transparent 
section,  and  2A  is  the  diameter  of  the  grating.  Since  in  the  manufacture 
of  circular  gratings  the  innermost  part  is  not  as  perfect  as  the  outer¬ 
most  parts,  the  first  four  periods  of  the  grating  have  been  ignored. 

This  experimental  necessity  will  turn  out  to  simplify  some  theoretical 


considerations . 


A  Qualitative  Explanation  of  Self-Imaging  with  Circular  Gratings 

Linear  gratings  are  known  to  self-image,  and  self-images  of 
circular  gratings  have  been  observed.  To  intuitively  understand  why 
circular  gratings  can  self-image,  consider  the  outermost  part  in  a  sector 
of  a  circular  grating.  In  this  region  the  circular  grating  is  approxi- 
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mately  like  a  linear  grating.  Therefore,  it  is  plausible  that  a 
circular  grating,  like  a  linear  grating,  will  self-image. 

Another  way  to  answer  this  is  to  show  that  the  circular  grating 
belongs  to  the  subset  of  radially  S5nranetric  objects  that  self-image. 

The  complete  set  of  self -imaging  objects  was  derived  by  Montgomery. 

He  assumed  that  a  plane  monochromatic  wave  impinged  upon  the  object.  He 
then  solved  the  wave  equation  in  Cartesian  coordinates  while  asking  the 
following  question:  **What  are  the  necessary  and  sufficient  conditions 
that  the  object  must  satisfy  in  order  that  a  faithful  image  of  it  be 
found  in  a  parallel  plane  z  =  Zrj.  >  0?”  The  solution  placed  restrictions 
on  the  spatial  frequency  domain  for  the  self-imaging  objects.  For 
objects  whose  spatial  frequency  is  much  less  than  (Montgomery  calls 

this  "weak  imaging")  we  have 

r=  2n/\z^\  =  /2n/Az^  (8) 

where  v  and  y  are  the  x  and  y  components  of  the  spatial  frequency,  n  is 
an  integer,  A  is  the  wavelength,  and  z-p  is  the  distance  for  self -imaging . 
Hence,  the  Fourier  spectrum  of  a  Talbot  object  can  only  exist  on  a  set 
of  concentric  circles  whose  radii  are  given  by  Therefore,  the  Fourier 

spectrum  for  the  subset  of  objects  with  radial  S3nnmetry  that  can  self- 
image  will  also  lie  on  these  concentric  circles.  By  changing  Montgomery's 
solution  to  polar  coordinates  and  using  the  fact  that  our  objects  have 
no  angular  dependence,  the  exact  solution  for  this  subset  can  be  deter¬ 
mined.  We  find  that  the  Fourier  spectrum,  in  addition  to  lying  on  the 
concentric  circles,  is  also  independent  of  the  angular  coordinate.  If 
the  circular  grating  belongs  to  the  Talbot  set  of  radially  symmetric 
objects,  then  its  Fourier  spectrum  must  behave  in  the  same  manner.  Since 
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Figure  8(a).  A  schematic  repre¬ 
sentation  of  a  circular  grating. 


Figure  8(b) .  The  transmittance  function 
of  the  circular  grating  where  the  first 
four  periods  have  been  ignored. 


Figure  9.  The  basic  Talbot  constant  radial  shearing  interferometer  setup. 


this  grating  can  be  expressed  as  a  Fourier  series  (eq.  (7)),  we  conclude 
that  its  Fourier  spectrum  is  probably  confined  to  equidistant  rings 
given  by  =  m/a.  These  rings  will  belong  also  to  the  larger  family  of 
non-equidistant  rings  of  eq.  (8)  if  we  select  n  =  0,  1,  4,  9...m  and 
a  =  Azt/2  or  =  2a^/X.  This  equation  is  the  same  as  the  original 
Talbot  formula  for  linear  gratings.  Hence  we  have  shown  qualitatively 
that  circular  gratings  are  a  subset  of  the  self-imaging  objects. 

Theory 

a.  General 

The  Talbot  interferometer  consists  of  two  identical  circular 
gratings  G1  and  G2  separated  by  a  distance  such  that  G2  lies  in  one  of 
the  negative  self-image  planes  of  G1  (see  Figure  9).  For  the  position 
of  G2  we  choose  a  negative  self-image,  midway  between  self-image  planes, 
so  that  diffevenoes  (and  not  sums)  in  the  object  will  be  observed.  Or 
we  obtain  these  differences  by  replacing  G2  with  a  negative  of  Gl.  The 
second  grating  is  now  placed  in  one  of  the  self-image  planes  of  the 
first.  We  also  use  a  plane  in  which  the  self-image  is  sharp  so  that  our 
fringes  will  have  their  highest  contrast.  The  gratings  are  aligned  with 
their  centers  coaxial  to  the  beam  direction.  Therefore  in  the  absence 
of  an  object  under  test  a  uniform  dark  field  is  observed  behind  G2. 
Between  the  two  gratings  at  a  distance  Z2  from  G2,  a  transparent  object 
is  inserted.  An  observer  focuses  his  eye  on  the  object  through  the 
second  grating  so  that  the  fringes  overlap  the  regions  to  which  they 
apply.  Although  a  white  light  source  may  be  used,  we  will  limit  our 
theoretical  considerations  to  a  monochromatic  plane  wave.  The  deriva- 
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tion  of  our  formulas  is  contained  in  the  next  section. 

A  monochromatic  plane  wave  impinges  upon  the  circular  grating  G1 
and  is  diffracted.  The  Huygens-Fresnel-Kirchhoff  equation^^  describes 
the  field  behind  G1  and  can  be  solved  by  the  method  of  stationary  phase^’^ 
The  solution  is  not  valid  in  the  inner  cone  of  diffraction  defined  by 
r  <  AzVa  (see  Fig.  10),  where  we  have  assumed  that  for  experimental 
reasons  the  first  four  periods  of  the  grating  are  ignored.  A  phenomeno¬ 
logical  justification  for  excluding  this  region  from  further  consideration 
in  the  instrument  lies  in  the  nature  of  the  Talbot  effect.  We  recognize 
that  outwardly  diffracted  rays  form  divergent  conical  wavefronts  and  in¬ 
wardly  diffracted  rays  form  convergent  conical  wavefronts.  The  inter¬ 
ference  between  a  plane  wave,  a  divergent  conical  wave  and  a  convergent 
conical  wave  is  utilized  for  the  circular  Talbot  effect.  These  three 
waves  overlap  in  the  shaded  areas  of  Fig.  11,  and  form  rings  around  the 
z-axis  with  a  ring  width  decreasing  to  zero  at  z  =  aA/2X.  The  inner  cone 
(cross-hatched  region)  is  one  of  the  zones  that  does  not  belong  to  the 
ring  domain.  This  region  is  to  be  ignored  in  evaluating  the  output  of 
the  instrument  and  is  usually  a  small  fraction  of  the  usuable  field. 

Ignoring  the  first  four  periods  of  the  grating  is  not  serious 
since  it  is  only  a  small  portion  of  the  total  grating  area.  We  could 
actually  block  off  this  region,  0  <  r  <  4a,  thereby  avoiding  contribu¬ 
tions  from  this  mathematically  unmanageable  area.  But  we  did  not  do  so 
since  we  know  that  light  from  that  small  portion  will  be  mostly  only  in 
the  inner  cone,  which  we  will  neglect  when  evaluating  the  output. 

The  result  of  our  analysis  in  the  next  section  gives  the  output 
of  the  instrument,  v(r,9),  as 
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Figure  10.  A  diagram  in  two  and  three  dimensions  illustrating  the  inner 
cone. 


G1 


Figure  11.  A  two  dimensional  drawing  of  the  propagation  of  the  first 

diffracted  orders.  The  Talbot  effect  occurs  only  within  the 
cross-hatched  regions. 
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C„  C„  _  („ 

(m)  (n) 

where  A  =  Xz^l^t  ^qM  is  the  object's  transmittance  function,  and 
^iTnn  ^  consequence  of  placing  G2  in  a  negative  Talbot  image.  This 
equation  contains  both  exponential  Fourier  series  of  G1  and  G2,  and 
shows  that  the  complex  amplitude  of  the  object  is  shifted  radially  by 
multiples  of  A.  The  amount  of  shear  depends  upon  the  distance  that  the 
object  is  placed  from  the  second  grating  (22)  as  well  as  the  angle  of 
the  diffracted  orders  (X/a). 

It  is  argued  in  the  next  section  that  further  simplification  of 
this  equation  can  be  achieved  by  observing  that  our  Fourier  coefficients 
decrease  by  1/m  and  that  terms  associated  with  high  frequencies  are 
decreased  by  whatever  observational  method  used.  Furthermore  C2>  C4,... 

=  0  if  we  use  square  grating  grooves.  The  result  of  these  approximations 
gives  us 

v(r,<p)  ~  Co2{uo(r,(p)  -  (Ci/Co)^  [uo(r  +  A,<p)  +  “  A.tp)])- 

As  with  the  case  of  linear  gratings  the  mark-to-space  ratio  is 
under  the  control  of  the  experimenter,  who  can  choose  (Ci/Cg)^  =  1/2 
(i.e.  mark-to-space  ratio  of  1.3),  so  that  v(r,  <p)  becomes 

v(r,9)  ~  uo(r,<p)  -  ■|■[uo(r  +  A,9)  +  uo(r  -  A,9)]. 

When  the  radial  shear  is  made  smaller  than  the  width  of  the  object 
details  this  equation  is  essentially  the  second  radial  derivative  of 
the  object, 

v(r,9)  ~  9^uo(r,<p)/3r^. 
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b«  Analysis 


The  Talbot  interferometer  with  circular  gratings  is  sketched  in 
Figure  9,  which  shows  all  the  appropriate  dimensions  and  symbols.  A 
plane  monochromatic  wave  impinges  on  Gl,  and  the  complex  amplitude  of 
the  wavefield  immediately  behind  Gl  is 


u(r,zi‘'')  =  exp(ikzi)  G(r), 


where  G(r)  is  the  transmission  function  of  the  circular  grating  given 

earlier.  We  will  ignore  all  constant  multiplicative  phase  factors  in 

this  derivation  since  the  final  observed  quantity  will  be  intensity. 

Beyond  Gl  we  have  Fresnel  diffraction  from  our  propagating  wave 

field  that  is  described  by  the  Huygens-Fresnel-Kirchhof f  integral, 

/•27T 


u(r,9,z)  T^l  I  uCr’jZi"*")  exp{i7r[r^  +  r’^  - 

2rr’cos(9  -  <P^)]/Az'}  r*dr*d9* 
where  (r,  9)  and  (r’,cp’)  are  polar  coordinates  in  plane  z  and  z^"^ 
respectively,  z’  =  z  -  zi,  and  the  approximation  sign  is  to  remind  us 
that  we  have  used  the  Fresnel  approximations  and  that  the  constant 
multiplicative  phase  factors  have  been  dropped. 

Integrating  over  cp',  substituting  G(r)  for  u(r,z'J’)  into  the 
remaining  integral  and  then  interchanging  summation  and  integration 
order  we  have 

u(r,z)  ~  ff. 

(m)  J 4a 

exp{iT7l2rm/a  +  (r*^  +  r^)/Az’]}  dr'. 

The  method  of  stationary  phase  that  solves  integrals  of  this 
type  cannot  be  applied  since  the  Bessel  function  oscillates  too  wildly 


-ml  Jo(27rrr '/Az') 
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over  the  region  of  statlonarity  to  be  approximated  by  a  constant.  By 
using  the  trigonometric  approximation  for  the  Bessel  function  we  can 
place  the  oscillating  parts  into  the  phase  term,  leaving  us  with  a 
slowly  varying  function.  The  approximation  is 

Jq(x)  ~  /2/i:x  cos(x  -  Tr/4);  x  >  25, 

and  therefore 

r'  Jo(2Trrr '/Az')  ~  /Az’r'/r  {exp  iTr[2rr'/Az'  -  1/4] 

+  exp  -iTr[2rf '/Az*  -  1/4]} 


if  2irrr*/Az'  >  25,  or  r  >  4Az'/r'. 

Recall  that  we  were  to  ignore  the  first  four  periods  of  the 
grating  for  experimental  convenience  since  the  innermost  parts  of  a 
circular  grating  are  more  difficult  to  manufacture  than  the  outer  parts. 
The  additional  value  to  this  exclusion  is  that  the  factor  Az'r '/r 
will  be  a  slowly  varying  function  for  r*  >  4a  and  that  therefore  the 
method  of  stationary  phase  is  applicable.  We  conclude  that  the  Bessel 
approximation  is  good  if  r  >  Xz'/a.  This  means  that  we  exclude  from  our 
consideration  in  the  instrument  a  region  defined  by  a  cone  about  the 


z-axis  whose  angle  with  the  axis  is  equal  to  the  angle  of  the  first 
diffracted  order  X/a.  Using  this  approximation  in  our  equation  for  the 
field  at  z  we  obtain 


u(r,z) 


;'r’/r  | 


exp  ITT 


(r'-r)^ 

Az' 


L  4-  * 

4  a 


+  exp  iir 


(r'+r)^  _  1  .  2mr ' 
Az'  4  a 


} 


dr’ 


This  equation  consists  of  two  parts;  the  first  part  accounts  for 
contributions  to  a  point  P  in  the  plane  z  from  the  grating  for  r  <0 
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giving  us  a  path  difference  of  (r  -  r*)  while  the  second  part  consists 
of  contributions  for  r'  <  0  that  give  a  path  difference  (r  +  r*).  Figure 
12  shows  these  two  contributions.  The  method  of  stationary  phase  is 
applied  to  the  integral. 

The  point  of  stationarity  for  the  first  part  is  ro'  =  r  -  mXz'/a. 
We  wish  to  know  the  complex  wave  field  at  P  that  lies  in  plane  z  at  a 
distance  r  from  the  optical  axis.  But  r  must  lie  within  the  regions 
where  self-imaging  interferences  occur  (see  Fig.  11).  Therefore,  not  all 
diffracted  orders  m  will  contribute  to  the  integral.  In  Figure  13  we 
have  plotted  the  stationary  point  as  it  falls  in  the  r  -  z  plane  for 
positive  values  of  m.  It  is  apparent  from  this  drawing  that  orders 
m  +  1  and  greater  will  fall  within  the  shaded  region  in  the  plane  z 
where  no  self-imaging  occurs.  Therefore  this  diagram  helps  us  to  deter¬ 
mine  the  positive  diffracted  orders  ^  m  that  place  ro*  within  the 
interval  of  integration  (4a,  A) .  In  like  manner  we  can  determine  the 
negative  diffracted  orders.  The  results  of  our  computation  give  us 

^  <  [(A  -  4a)a/X](l  +  Imj)  =  (N  -  4)z-p/2(l  +  Imj);  A  =  Na 

where  zj  =  2a^/X  is  the  self-imaging  distance  and  N  is  the  number  of 
lines  in  the  grating. 

When  z  is  a  small  fraction  of  (N/2  -  2)zrj.  then  many  orders  of  m 
contribute.  The  strength  of  the  various  orders  is  proportional  to  the 
Fourier  coefficients.  As  we  have  stated  earlier,  only  the  zeroth  and  two 
first  orders  are  significant.  In  this  case  z  <  Nzx/4  -  zj  ^  Aa/2X,  as 
stated  earlier.  We  will  carry  all  the  orders  in  our  computation  and  take 
only  those  into  account  that  lie  within  the  instrument’s  operating  range 
when  evaluating  the  output. 
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Figure  13.  A  plot  of  the  stationarity  point  for  that  region  of  the 
grating  that  lies  above  the  x— axis  (r  —  r'). 
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Applying  the  method  of  stationary  phase  to  the  first  part  gives 
us 

u,(r,:,)  »  A  -  (mX^'/ra) 

(m) 

where  ui(r,z)  is  the  contribution  due  to  the  first  part  of  the  integral. 

The  point  of  stationarity  for  the  second  part  of  the  integral  is 
to’  =  -  r  -  mAz'/a.  For  m  ^  0,  the  stationary  point  lies  outside  of  the 
region  of  integration  (4a,  A),  and  therefore  does  not  contribute  to  the 
integral.  For  m  =  -1,  the  stationarity  lies  inside  the  inner  cone  that 
we  are  excluding  from  our  instrument  in  order  that  our  approximation  be 
valid.  The  remaining  values  of  m  (m  <  -1)  do  not  significantly  contri¬ 
bute  to  the  integral  since  this  contribution  is  multiplied  by  the  associ¬ 
ated  Fourier  coefficients  which  are  decreasing  by  1/m,  or  by  (1/m)^  in 
intensity.  Thus  the  contribution  of  this  second  part  to  the  integral  is 
negligible  in  comparison  with  the  first  part,  and  therefore  our  field 
u(r,z)  is  given  by  the  computation  for  ui(r,  z) 

u(r,z)  «  ui(r,z) 

«  ^  Cj,  /l  -  (mAzVra)  ^  ^2Trimr/a^ 

(m) 

In  the  plane  just  preceding  the  object  (z  =  0“) ,  the  field  is 

u(r,0-)  ~  A"+"  (mAzi/ra)  ^  ^2Trimr/a^ 

(m) 

and  just  behind  the  object  we  have 

u(r,cp,0''‘)  =  u(r,0“)  uo(r,«p), 

where  uo(r,<p)  is  the  complex  amplitude  of  our  object  that  we  want  to 
test.  In  the  region  beyond  the  lens  the  propagating  wave  field  is 
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computed  again  by  use  of  the  Huygen— Fresnel— Kirchhoff  formulation, 

Jl^fy(z)  r2Tr 

I  u(r,^,0'^) 

f^(z)  Jo 

exp{i7T[r^  +  -  2rr*cos(<p  -  <P^)]/^z} 

r ’ dr ’ d9  • , 

where  we  limit  the  integration  to  the  self-image  field,  f^Cz)  =  A  -  Az/a, 
and  f£(z)  =  4a  +  Xz/a.  We  expand  our  test  object  uo(r,9)  in  a  Fourier 
exponential  series 

uo(r,9)  = 

(P) 

since  any  object  is  periodic  in  9.  Placing  this  into  our  integral  and 
integrating  over  9',  we  have 

iTT(m/a)^Xzi  ip(9-Tr/2) 


u(r,9,z) 


If  EE 

(m)  (p) 


Cm  ^ 


‘fu(z) 


'f»(z) 


r'/l  +  mXzi/r'a  Jp(27rrr'/Xz)  u  (r’) 


exp{iTr[(r^  +  r'^)/Xz  +  (2mr'/a)]} 


In  order  to  solve  this  by  the  method  of  stationary  phase  we 
approximate  the  Bessel  function  as  before  and  again  find  that  we  must 
ignore  the  inner  cone.  Since  this  new  inner  cone  lies  within  the  inner 
cone  of  the  grating  Gl,  no  further  restrictions  are  required  in  the 
theory  of  the  instrument.  Hence, 


2Trr'  /l  +  mXzi/r’a  J„(2Tirr '/Az)  ~ 


/Xz/r  v^r"'  +  mXzi/a  cos[2Trrr*/Xz  -  Tr/4  —  pTr/2], 
where  r  >  Xz/a  and  r'  >  4a  as  before. 

Let  the  cosine  term  be  represented  by  its  exponential  equivalence. 
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and  our  field  equation  becomes 


u(r,9,z)  ~  A/Azr 


_  iir(m/a)2Xzi  ip((p-Tr/2) 
e  e 


(m)  (p) 


’fu(z) 


'fo(z) 


'x'  +  mXzi/a 


{exp  iTT[(r'  -  r)^/Az  +  1/4  +  p/2  +  2mr'/a]  + 
exp  iTr[(r'  +  r)^/Az  -  1/4  -  p/2  +  2mr'/a]}  dr' 

The  second  part  is  approximately  zero  as  before,  and  the  first  part 


gives  us 


u(r,cp,z) 


EE 

(m)  (p) 


-iTr(m/a)2A(z-Zi  )  r - - x-, - 

Cjn  e  ^  /I  -  mA(z  -  z^j/ra 


u.(r  -  mxz/a)  .‘W" 


Note  that  our  test  object  is  duplicated  and  shifted  by  multiples 
of  Xz/a.  The  field  just  prior  to  the  second  grating  is  obtained  by 
setting  z  =  22*  The  above  equation  becomes 


u(r,cp,Z2“) 


~E 


-  Zi)/ra 


uo(r  -  mAz2/a,9)  e 


2Trimr/a 


where  we  have  used  uo(r,(p)  =  7  u_(r)  e 


The  field  just  behind  the  second  grating  is 


u(r,9,Z2'^)  =  u(r,(p,Z2  )  G(r)  «  -  mX(z2~zi)/ra 

(m)  (n) 

“i7rX(m/a)^(z2-zi)  .  a  /  27Ti(m+n)r/a 

e  uo(r  -  mXz2/a,9)  e  .  .  (10) 


For  maximum  contrast,  we  choose  Z2  -  zi  =  Ma^/X,  where  M  is  an  integer. 
Furthermore,  since  we  want  to  see  diffevenoes  (not  sums)  in  object 
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detail,  we  choose  that  M  should  be  an  odd  integer.  We  have  exp  iTTX(m/a)^ 
(z^  -  zi)  =  exp  iirm,  which  gives  us  the  necessary  negative  sign  for  odd 
integers.  The  full  significance  of  this  is  deferred  until  the  spatial 
filtering  section.  With  this  distance  our  equation  becomes 


u(r,<p,Z2‘*')  ~ 

Cji  e^™  ^2iri(m+n)r/a  mMa/r  uo(r  -  mXz2/a,<p). 

(m)  (n) 

This  equation  contains  both  exponential  Fourier  series  of  G1  and 
G2,  and  shows  that  the  complex  amplitude  of  the  field  is  shifted  radially 
by  multiples  of  Xz2/a.  The  amount  of  shear  depends  upon  the  distance 
that  the  object  is  placed  from  the  second  grating  as  well  as  the  angle 
of  the  diffracted  orders.  The  square  root  term  accounts  for  the  increase 
or  decrease  in  the  size  of  the  object  due  to  either. 

To  simplify  our  mathematical  manipulations  we  would  like  to  work 
in  the  area  where  /T  ±  Ma/ r  ~  1  ±  Ma/2r  ^1.  This  means  that  Ma/2r  <<  1 
or  r  »  Ma/2.  Recall  from  eq.  (10)  that  this  term  is  /l  -  X(z2  -  zi)/ra. 
Therefore,  r  >>  X(z2  -  zi)/2a.  This  states  that  r  must  be  greater  than 
a  cone  which  lies  within  the  inner  cone  that  we  are  already  ignoring  for 
aforementioned  reasons.  Thus  the  square  root  term  is  close  to  1 
for  the  region  outside  the  inner  cone.  Thus  our  equation  slmpli~ 

fies  to 

v(r,cp)  =  u(r,<p,Z2+) 


-EE 

(m)  (n) 


Cm  Cl 


imn  2TTi(m+n)r/a 


uo(r  -  mA,(p), 


(11) 


where  A  =  Xz2/a  is  the  shear  introduced,  and  v(r,q))  is  a  new  s)anbol  to 
indicate  the  observed  field. 
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The  only  significant  terms  of  this  equation  are  for  m  = 

Our  Fourier  coefficients  decrease  by  a  factor  1/m.  And  the  next  signi¬ 
ficant  coefficient  C3  (€2^  0)  is  one-third  that  of  .  Since  our 
Fourier  coefficients  are  in  pairs,  we  can  compare  C3C0,  which  is  the 
next  largest  term  that  we  can  ignore,  with  which  is  the  smallest 

term  that  we  will  keep,  and  see  that  the  term  that  we  will  throw  away  is 
less  than  one-half  that  of  the  smallest  term  we  keep.  Furthermore  this 
n  =  0,  m  =  3  term  contains  fairly  high  spatial  frequencies  as  evidenced 
by  the  term  exp  27ri3r/a.  These  high  frequencies  will  be  further  reduced 
in  contrast  by  whatever  observational  method  is  used:  eye,  camera,  or 
film.  A  detailed  justification  of  this  statement  is  involved  since  we 
would  have  to  consider  the  spatial  frequency  spectrum  of  the  intensity 
(after  a  modulus  square  operation).  Under  these  assumptions  our  equation 
reduces  to 

v(r,(p)  ~  uo(r,<p)  -  Ci  2{uo(r+A,(fi)  +uo(r-A,(p)} 

+  CgCi  cos(2Trr/a)  uo(r,9) 

-  Uo(r+A,9)  e  +  CgCi  e  \ 

-  uo(r-A,9)  e  +  CgCi  e  } 

As  we  have  just  stated,  the  terms  with  the  exponential  quantities  contain 
very  high  frequencies  and  the  contrast  of  these  terms  will  be  attenuated 
or  even  eliminated  by  the  observational  method  as  with  a  grating  with  a 
period  too  fine  for  the  eye  to  resolve.  Therefore, 

v(r,9)  ~  Cq^  uo(r,cp)  -  Ci^  uo(r+A,<p)  +  uo(r-A,(p)  . 

If  (Ci/Cq)^  =  1/2,  then  we  have 

v(r,cp)  ~  Co^{uo(r,<p)  -  l/2[uo(r+A,9)  +  uq  (r-A,cp)  ]  } , 
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and  this  equation  is  approximately  the  second  radial  derivative  when  the 
shear  is  smaller  than  the  object  details  of  interest.  Hence 

v(r,<p)  ~  92uo(r,(p)/3r^. 

To  obtain  (Ci/Cq)^  =  sinc^a  =  1/2  we  choose  a  =  0.44  or  a  mark- 
to-space  ratio  of  1,3.  In  conclusion,  we  see  that  the  output  of  the 
Talbot  interferometer  with  circular  gratings  is  triple  shearing  inter¬ 
ferences.  These  triple  shearing  interferences  become  the  second  deriva¬ 
tive  when  the  shear  is  made  smaller  than  the  object  details.  The  shear 
in  turn  is  proportional  to  the  distance  that  the  object  under  test  is 
placed  from  the  second  grating. 


Shifted  Circular  Gratings 

In  this  technique  the  second  grating  which  lies  in  one  of  the 

self-image  planes  is  shifted  laterally  a  distance  s  from  the  self-image. 

Thus  the  rings  of  equal  diameter  for  G2  and  the  self-image  overlap  on 

the  perpendicular  bisector  of  the  line  joining  their  centers.  A  phase 

object  inserted  between  the  gratings  a  distance  z  from  G2  causes  the 

rays  passing  through  it  to  bend  by  an  angle  ot  (see  Figure  14(a)).  Thus 

the  rings  of  the  self-image  are  distorted  and  the  overlapping  rings  will 

deviate  from  the  perpendicular  bisector  as  in  Figure  14(b).  We  will 

18 

generalize  Lau’s  Dupligramm-method  from  Fresnel  Zone  Plate  rings  to 
equidistant  rings.  The  deviation  will  be  shown  to  measure  the  aberration 
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Figure  14(a),  A  ray  diagram  Figure  14(b).  The  two  circular  gratings  are 
illustrating  how  an  object  shifted  with  respect  to  each  other  causing 

bends  the  rays  in  a  circular  moire  fringes.  When  one  circular  grating  is 

Talbot  interferometer.  deformed  then  the  distance  from  the  y-axis 

of  the  intersections  is  a  measure  of  this 
distortion. 


Figure  15.  Using  the  lens  under  test  in  collimation  increases  the  sensi¬ 
tivity  of  the  test  since  the  unwanted  focusing  term  is  biased 
out  in  collimating  the  source. 


of  a  test  lens,  where  a  =  a(r)  describes  the  angular  aberrations  of  the 
lens  under  test. 

The  opaque  circular  rings  of  G2  lie  on  radii  =  ma  (a  is  the 
grating  period).  The  distorted  rings  of  the  self-image  lie  on  radii 
rjjj’  =  rjj^  -  az.  We  choose  our  coordinate  axis  in  the  plane  of  G2  such 
that  the  x-axis  joins  the  two  centers  of  the  two  concentric  circular 
ring  structures  and  the  y-axis  lies  on  the  perpendicular  bisector  as 
shown  in  Figure  14(b).  The  equation  for  the  circles  of  our  two  ring 
structures  are: 

(x  +  s/2)2  +  y2  =  (x  -  s/2)2  +  y2  = 

The  simultaneous  solution  of  these  two  equations  gives  us  the 
coordinates  of  the  overlapping  rings.  Solving  for  6(rin),  which  is  the 
x-coordinate  of  the  rnoir^  fringe  closest  to  the  y-axis,we  have 

6(rj^)  =  [a(%)z/s]  [rm  -  a(rm)z/2]  . 

We  can  always  choose  the  lateral  shift  a(rin)z  <<  2rj„  so  that  the  devia¬ 
tion  is  simply  related  to  the  deviations, 

5(rm)  «  z  rm  a(rin)/s. 

As  an  example  we  pick  a  lens  where  +  ...  .  The  term 

a^r  describes  the  focal  power  while  a3r^  +  ...  are  the  spherical  aberra¬ 
tion  terms  of  interest.  We  can  eliminate  the  air  term  by  using  the  test 
lens  in  collimation  as  shown  in  Figure  15.  In  this  configuration  we  have 
a,  -  air  =  a(r„).  Therefore 

6(rm)  =  [a(rj„)z/s]  [rj^  -  a(rni)z/2]. 

For  small  distances  r„,  a(rn,)  «  0  since  the  lens  has  no  aberrations  on 
axis.  For  larger  distances  r^j  we  choose  a(rjj)z  <<  2rjjj  so  that 
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^  (zr^/s)  ot(rjjj). 

The  distance  z  is  equal  to  the  separation  of  the  gratings.  We  see  that 
the  deviations  from  the  y-axis  oi(r^)  are  simply  related  to  the  angular 
aberrations  of  the  lens. 


Filtering 

The  setup  for  spatial  filtering  is  shown  in  Figure  16.  It 
requires  the  addition  of  a  telecentric  system  placed  behind  G2  to  make 
it  possible  to  separate  out  the  various  diffracted  orders,  by  placing  a 
binary  mask  in  the  filter  plane.  The  analysis  and  setup  are  analogous 
to  the  lateral  shearing  interferometer.  The  results  are  identical  if  we 
replace  the  variable  x  by  r.  The  usual  bandwidth  restrictions  for 
spatial  filtering  apply.  Thus  the  bandwidth  Ap  of  our  object  UQ(r,  cp) 
must  be  less  than  1/a.  Physically  this  means  that  the  grating  period 
must  be  less  than  the  tiniest  object  details  that  we  wish  to  resolve. 

a.  Zeroth-order  Filtering 

By  zeroth-order  filtering  we  mean  that  only  those  rays  that  are 
approximately  parallel  to  the  optical  axis  beyond  G2  are  allowed  to 
reach  the  observation  screen.  This  is  accomplished  with  a  telecentric 
system  that  focuses  these  rays  to  a  point  on  the  optical  axis.  This  is 
illustrated  in  Figure  17.  A  binary  filter  with  an  aperture  large  enough 
to  allow  these  rays  to  pass,  but  small  enough  to  reject  the  other  orders, 
is  placed  in  this  plane.  Thus  only  terms  with  m  +  n  =  0  in  eq.  (10) 
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Figure  16.  The  Talbot  interferometer  modified  to  do  spatial  filtering. 

The  telecentric  system  provides  a  filter  plane  for  the 
insertion  of  the  binary  filter  and  the  image  plane  is  con¬ 
jugate  to  the  object. 


L2 


Figure  17 . 


A  diagram  showing  how  the  zeroth  diffracted  order  is  filtered. 
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will  appear  in  the  final  output  of  the  instrument.  In  order  that 
filtering  be  accomplished,  the  zeroth  order  must  not  overlap  terms  with 
different  (m  +  n)  orders.  Thus  the  bandwidth  Ap  of  our  object  UQ(r,  (p) 
must  be  less  than  the  inverse  of  the  period  of  our  grating  (Ap  <  1/a). 
Physically,  this  means  that  the  grating  period  must  be  smaller  than  the 
tiniest  object  details  that  we  wish  to  resolve.  If  we  place  an  opaque 
mask  with  a  disc  opening  of  radius  Af/2a  into  the  back  focal  plane  of  L2, 
then  what  we  observe  at  the  screen  is 

v(r,9)  ~  ^  uo(r  -  mA,(jp). 

(m) 

Only  terms  with  m  =  ~1, 0,4-1  of  this  equation  are  significant,  since  C2 
^  0  and  ^  Ci^/9.  Therefore,  we  have 

v(r,<p)  ~  Uo(r,9)  -  (Ci/Cq)^  [ugCr+A.cp)  +  uo(r-A,9)]  . 

If  we  choose  a  mark-to-space  ratio  of  1.3,  then  (Ci/Cq)^  =  1/2, 
and  the  above  equation  becomes 

v(r,(p)  ^  Uo(r,(p)  -  (1/2)  [uo(r4A,<p)  +  uq (r-A,(p) ] . 

If  the  shift  A  is  much  less  than  the  tiniest  object  details  we  wish  to 
observe,  then  this  equation  approximates  the  second  derivative,  v(r,cp) 
3^UQ(r,(p)/Br^.  This  result  is  similar  to  our  previous  unfiltered 
result.  There  we  did  low-pass  filtering  with  the  frequency  limitations 
of  the  observational  method.  When  the  gratings  can  be  resolved  by  the 
eye  or  other  recording  media,  then  one  sees  the  grating  structure  super¬ 
imposed  on  the  object.  For  this  case,  zeroth  order  filtering  would  be 
desirable. 
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b.  First-order  Filtering 


First-order  filtering  is  just  as  easily  accomplished,  but  the 
setup  must  be  modified  slightly  in  order  to  achieve  the  first  radial 
derivative.  Figure  18  demonstrates  that  both  first-order  diffraction 
rays  intermingle  in  the  back  focal  plane  of  the  lens,  making  it  impossi¬ 
ble  to  separate  out  either  the  m  +  n  =  +1  or  -1  term.  A  simple  remedy 
to  restore  order  is  to  block  one  half  of  the  grating.  In  this  way  rays 
from  the  lower  half  do  not  intermingle  with  the  upper  rays  in  the 
diffraction  pattern. 

Since  the  m  +  n  =  +1  or  -1  term  lies  in  a  semicircle,  an  opaque 
mask  with  a  transparent  semi-annulus  will  transmit  either  one  of  the 
first  orders  (see  Fig.  19).  This  modified  setup  means  that  one  half  of 
our  circular  grating  is  blocked. 

Fortunately,  our  theoretical  analysis  is  only  modified  by  a 
constant  factor.  1/2  as  discussed  in  Appendix  3.  Thus  the  output  of 
this  modified  setup  is  1/2  that  of  eq.  (9),  i.e.. 


v(r,9) 


(m)  (n) 


(m+n)r/ a 


uo(r  -  mA,<p). 


With  the  insertion  of  the  semi-annular  binary  mask,  only  the 
m  +  n  =  +1  (or  m  +  n  =  -1,  if  we  rotate  the  annulus  by  tt  radians)  will 
pass  through  to  the  screen.  With  a  Ronchi  ruling  all  even  Fourier 
coefficients  are  zero.  Therefore,  the  field  at  the  screen  is 

v(r,cp)  ~  yCoCi  [uo(r,9)  -uo(r-A,9)] 


If  we  drop  the  uninteresting  factor  exp  2irir/a  from  further 
consideration,  and  Introduce  a  shear  that  is  smaller  than  the  tiniest 
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Figure  18.  A  diagram  showing  the  intermingling  of  the  negative  and  positive 
first  diffracted  orders,  thereby  making  it  impossible  to  do 
first  order  filtering. 


Figure  19.  The  tail  end  of  the  interferometric  setup  showing  the  semiannular 
spatial  filter  used  with  first  order  filtering. 
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object  details  of  interest,  then  our  equation  reduces  to  the  first 
radial  derivative,  v(r,9)  ^  8u(r,9)/9r. 


Experiments 

Our  circular  gratings  were  photographic  reductions  of  the  circu¬ 
lar  pattern  that  is  included  in  the  Edmund  Scientific  moire  kit.  Two 
identical  gratings  were  made.  Our  gratings  were  examined  under  a  micro¬ 
scope  and  determined  to  have  a  mark-to-space  ratio  of  about  1 . 

Spatial  filters  were  made  from  photographs  of  drawings  that  were  reduced 
to  the  required  size.  Both  xenon- arc  and  helium— neon  laser  were  used  as 
sources.  Color  fringes  of  high  contrast  were  observed  with  the  xenon 
arc.  We  present  only  our  black  and  white  work. 

The  interferometer  was  set  up  on  a  laboratory  bench  according  to 
the  schematic  of  Figure  9.  A  5  mW  helium-neon  laser  beam  was  expanded 
to  50  mm  by  a  250  mm  effective  focal  length  objective  and  filtered  with 
a  6-micron  pinhole.  Our  first  circular  grating  G1  was  placed  within  the 
beam  and  another  identical  grating  G2  was  centered  in  one  of  the  negative 
self-image  planes  of  Gl.  The  negative  self-image  plane  in  which  G2  is 
inserted  is  chosen  so  that  the  test  object  which  is  between  the  gratings 
can  be  placed  sufficiently  far  away  from  G2  to  yield  the  desired  shear. 

A  uniform  dark  field  is  observed  at  the  screen.  Inserting  a  test  object 
produces  moire  fringes  through  a  distortion  of  the  first  grating’s 
negative  self-image  at  G2.  For  photographic  purposes,  a  lens  was  placed 
behind  G2  with  the  film  plane  conjugate  to  the  object. 
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a.  No  Filtering 

The  fringe  pattern  for  a  bifocal  lens  is  shown  in  Paper  2,  Fig. 

2.  We  used  Kodak’s  high-contrast  copy  film.  The  fringes  are  contours 
of  constant  phase  derivative.  Note  that  the  smaller  lens  contains  more 
fringes  than  the  larger  one  due  to  its  higher  focal  power.  The  object  is 
aligned  so  that  the  optical  axis  for  the  large  lens  is  collinear  with  the 
optical  axis  of  the  system  as  we  see  from  the  circular  fringes.  The 
smaller  lens  has  horseshoe-type  fringes  since  it  is  off-axis  and  we  are 
looking  at  its  off-axis  aberrations. 

Moire  fringes  in  the  shape  of  hyperbolas  are  observed  when  two 
circular  gratings  are  shifted  laterally  with  respect  to  each  other,  and 
look  very  much  like  spikes  for  small  shifts.  If  our  test  object  is  a 
candle  flame,  the  spikes  will  be  deformed  by  the  change  in  refraction 
index  of  air  due  to  heat.  In  Paper  2,  Fig.  3  is  a  photograph  of  this 
experiment  where  the  candle  is  approximately  13  cm  from  G2. 

b.  Shifted  Gratings 

As  stated  above,  shifted  circular  gratings  produce  moire  fringes 
in  the  shape  of  hyperbolas.  There  is  also  a  moire  fringe  due  to  the 
overlapping  of  rings  of  equal  diameter  that  fall  on  the  perpendicular 
bisector  of  the  line  joining  the  centers  of  the  shifted  gratings.  We 
can  display  this  fringe  alone  by  using  a  small  shift  s  between  the 
grating  centers. 

It  was  shown  earlier  that  with  the  addition  of  a  phase  object 
the  deviation  6  from  the  perpendicular  bisector  is  simply  related  to  the 
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angle  through  which  the  ray  is  bent.  For  the  case  of  a  lens  in  collima- 
tion,  6  was  simply  related  to  the  aberrations  of  the  lens.  Figure  20 
shows  a  lens  of  380  mm  focal  length  in  the  collimation  setup.  The  shift 
s  =0.07  mm  and  z  =  57  cm.  The  deviations  ^  are  a  measure  of  this  lens’s 
aberration. 

c.  Filtering 

Filtering  is  accomplished  by  means  of  binary  masks  placed  in  the 
Fourier  plane  of  the  transform  lens  L2  that  separates  one  diffracted 
order  from  another.  The  location  of  the  filter  plane  depends  on  the 
focal  length  of  L2  and  sometimes  also  on  the  object  if  it  happens  to  be 
a  lens.  Should  the  test  object  be  a  bifocal  lens  with  focal  lengths  f^ 
and  f25  for  example,  then  two  Fourier  planes  exist  and  filtering  cannot 
be  done  in  either  plane  without  eliminating  information  located  in  the 
other.  Under  certain  conditions  there  exists  an  intemediate  plane  in 
which  the  orders  are  separated,  and  the  filter  can  be  inserted  there. 

The  condition  under  which  this  region  exists  is  given  in  Appendix  4  as 
A/a  >  where  h^  is  the  height  of  that  part  of  the 

object  with  focal  length  f^,  y  is  a  number  less  than  1  and  depends 
upon  the  angle  at  which  the  diffracted  orders  meet,  and  M  is  the  magni¬ 
fication  of  the  object  at  the  screen.  Note  that  the  condition  is 
independent  of  the  focal  length  of  the  transforming  lens. 

Our  demonstration  of  triple  and  double-shearing  interferences 
for  zeroth  and  first-order  filtering  respectively  uses  a  ring  cut  from 
brass  tubing  with  an  outside  diameter  of  19.05  mm  and  an  inside  diameter 
of  15.37  mm  and  soldered  to  a  post.  If  the  shear  A  exceeds  the  thickness. 
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then  individual  rings  will  be  observed.  In  this  way  it  is  sufficient  to 
count  the  number  of  rings  to  determine  whether  double  or  triple  shearing 
interferences  are  obtained  by  simple  filtering.  The  ring  was  inserted  a 
distance  973  mm  from  G2.  That  gave  a  shear  (A  =  Az2/a.)  of  2.4  mm,  there¬ 
by  insuring  that  our  images  would  be  separated. 

For  zeroth-order  filtering  a  binary  mask  with  a  transparent  disc 
is  placed  in  the  filter  plane.  The  diameter  of  the  disc  is  such  that 
the  periphery  of  the  disc  falls  midway  between  the  zeroth  and  first 
diffraction  orders.  At  the  observation  screen  we  observe 

|v(r,cp)|2  ^  Cq**  |uo(r,cp)  I  2  +  luo(r+A,9)  p  +  .|uo(r-A,(p)  |2}. 

The  intensity  of  the  unshifted  ring  is  associated  with  Cq^  while  the 
intensity  of  the  shifted  rings  are  associated  with  Ci^.  Hence  the 
shifted  rings  will  be  of  lower  contrast  than  the  unshifted  ring.  A 
photograph  of  this  is  presented  in  Fig.  21(a).  One  half  of  the  field  is 
blocked  for  convenience  when  comparing  this  figure  with  the  first-order 
filtering  of  Fig.  21(b). 

For  first-order  filtering  one  half  of  the  circular  grating  is 
blocked  to  separate  positive  and  negative  orders;  a  semiannular  ring  is 
placed  in  the  filter  plane,  oriented  to  pass  either  n+m=+lorn+m 
=  -1  diffracted  orders.  In  the  photograph  of  Figure  21(b),  we  show  the 
result  of  passing  the  positive  order.  Two  rings  of  equal  contrast  are 
observed  as  predicted  by  the  theory.  One  image  is  unshifted  and  the 
other  is  shifted  away  from  the  axis.  Rotating  the  semiannulus  by  180 
degrees  allows  the  negative  diffracted  order  to  pass.  In  this  case,  we 
would  have  an  unshifted  image  and  another  image  of  equal  intensity 
shifted  toward  the  axis.  Comparing  this  photograph  with  that  of  the 
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previous  figure  demonstrates  that  spatial  filtering  with  simple  binary 
masks  can  produce  either  double  or  triple  shearing  interferences. 


Figure  20.  The  aberrations  of  the  lens  tested  with  shifted  circular 
gratings  are  directly  displayed  as  deviations  from  a 
straight  line  passing  through  the  center.  Virtually  no 
aberrations  are  visible  near  the  optical  axis. 
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Figure  21(a) .  A  ring  object  under  test  in  the  constant  radial  interfero¬ 
meter.  A  large  shear  is  introduced  by  placing  the  object 
93  cm  from  G2.  Zeroth  order  filtering  produces  triple 
shearing  interferences  and  only  one  half  of  the  field  is 
viewed . 


Figure  21(b) .  The  same  setup  as  (a)  except  that  first  order  filtering  is 
used . 
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CHAPTER  IV 


THE  TALBOT  SPECTROMETRIC  IMAGER 

Introduction 

Of  the  many  instruments  that  measure  the  spectrum,  Fourier  spec¬ 
troscopes  have  received  considerable  attention.  This  is  due  mainly  to 
two  advantages:  the  multiplex  advantage  and  the  optical  acceptance 
advantage.  The  latter  applies  to  the  Michel son  interferometer  only  and 
theoretical  luminosity  gains  of  200  over  the  grating  spectrometer  have 
been  reported.  The  multiplex  advantage  refers  to  the  simultaneous 
encoding  of  all  wavelengths  onto  a  single  detector  by  assigning  each 
resolution  element  its  own  carrier  frequency.  In  this  manner  a  time 
advantage  of  great  importance  is  obtained  which  is  used  either  to  decrease 
the  time  of  observation  or  to  increase  the  signal— to— noise  ratio  by  time 
averaging.  An  example  of  this  time  advantage  is  the  planetary  spectra 
of  P*  Connes  and  J.  Connes  obtained  with  a  Michelson  interferometer  that 

could  have  been  obtained  with  a  good  grating  spectrometer  in  a  total 

2 

time  of  a  few  thousand  years. 

In  addition  to  the  Michelson  interferometer,  various  additional 
schemes  for  multiplexing  a  spectrum  have  been  devised.  The  output  of  a 
dispersing  spectroscope  such  as  a  prism  has  been  multiplexed  in  two  ways. 
One  method  uses  a  rotating  reticle  on  which  a  set  of  concentric  annular 
rings  are  marked  and  places  it  in  the  exit  plane  of  the  spectroscope.^ 

Each  annular  ring  has  a  different  number  of  equal  opaque  and  transparent 
sections  and  encodes  a  different  spectral  element.  The  total  light  is 
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gathcrsd  and  recordad  •  Tlia  othar  inGthod  is  to  placG  a  stationary  mask. 

4 

in  the  exit  plane  of  the  dispersing  spectroscope  .  The  mask  has  a  number 
of  slits  that  are  either  opaque  or  transparent.  The  light  is  gathered  and 
recorded,  and  then  a  new  mask  is  inserted.  The  number  of  resolvable 
spectral  elements  equals  the  number  of  masks  that  are  used.  In  order  to 
retrieve  the  spectrum,  the  masks  are  designed  so  that  the  recorded  inten- 
sities  can  be  expressed  as  a  product  of  a  Hadamard  matrix  and  the  unknown 
spectral  elements.  Recovery  is  accomplished  by  a  Hadamard  transformation 
which  is  reportedly  five  times  faster  on  the  computer  than  the  Fast 
Fourier  transform.  Of  importance  in  this  chapter  is  the  scheme  for 
multiplexing  spectra  with  the  Talbot  effect."^  The  coding  is  a  cosine 
transformation  of  the  spectrum  as  a  function  of  the  separation  of  two 
gratings.  Another  cosine  transformation  recovers  the  spectrum. 

Not  only  can  the  spectrum  be  multiplexed  but  images,  too.  An 
image  scanning  system  that  multiplexes  a  two  dimensional  scene  with  a 
single  detector  has  been  proposed.  The  scene  is  imaged  onto  a  two 
dimensional  mask  that  either  transmits  or  blocks  light  at  various  points 
on  its  surface.  The  light  is  gathered  and  the  intensity  is  recorded. 

Then  this  mask  is  replaced  by  another  mask  with  different  opaque  and 
transparent  regions  than  the  previous  one,  and  the  intensity  is  recorded 
again.  If  m  x  n  elements  of  the  scene  are  to  be  encoded  then  m  x  n  masks 
QX&  prepared  and  m  x  n  intensities  recorded.  The  masks  are  based  on  the 
two  dimensional  Hadamard  matrices  and  retrieval  is  accomplished  with  a 
Hadamard  inverse  transformation. 

This  image  scanner  was  extended  to  encode  the  brightness  distri¬ 
bution  of  a  scene  not  only  as  a  function  of  coordinates  (a, 3)  as  above 
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but  also  as  a  function  of  A  from  point  to  point.  This  new  instrument  is 
called  a  spectrometric  imager.  Thus  the  decoded  output  of  such  a  system 
could  be  a  series  of  scenes  in  which  the  intensity  at  a  given  A  would  be 
presented.  The  modification  here  is  to  follow  the  mask  of  the  above 
image  scanning  system  with  a  dispersing  spectrometer  like  a  prism  and 
then  a  slotted  mask  with  say  p  slits  placed  in  the  exit  plane  of  it.  The 
slits  would  be  designed  according  to  the  Hadamard  matrix  again;  the 
light  gathered;  and  the  intensity  recorded.  Here  m  x  n  x  p  elements  are 
encoded  which  require  m  x  n  x  p  measurements. 

The  Talbot  interferometer  follows  this  kind  of  development. 

First  there  was  the  spectrometer  and  now  in  this  chapter  we  examine  the 
theory  of  using  the  effect  to  make  an  image  scanner  and  a  spectrometric 
imager.  Our  theory  for  the  image  scanner  shows  that  the  spatial  distri¬ 
bution  of  a  scene  is  encoded  as  a  Fourier  cosine  transform  if  we  restrict 
the  acceptance  angle  of  the  instrument.  The  scene  is  recovered  from  the 
encoded  data  by  a  cosine  transformation.  The  Talbot  equivalent  of  a 
spectrometric  imager  is  a  bit  more  complicated,  requiring  four  gratings 
if  the  total  Fourier  cosine  space  is  to  be  accessible.  The  two  grating 
instrument  can  only  explore  a  conical  surface  in  the  Fourier  space  of 
S(a,g,A).  The  extra  set  of  gratings  allows  us  to  explore  the  whole 
volume.  The  theory  for  the  Talbot  image  scanner  and  the  Talbot  spectro¬ 
metric  imager  is  presented  in  this  chapter.  The  image  scanner  is  a 
special  case  of  the  spectrometric  imager. 
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Theory 


The  basic  theory  for  an  image  scanner  and  a  spectrometric  imager 
(See  Fig.  22)  consists  of  two  gratings  and  resembles  the  shearing  inter¬ 
ferometers  of  Chapter  III.  It  differs  in  the  way  we  observe  the  light 
intensity.  The  light  behind  the  second  grating  is  collected  onto  a 
single  detector,  and  its  intensity  as  a  function  of  grating  separation  is 
measured.  The  second  grating  has  a  lateral  motion  that  allows  electronic 
filtering  of  the  encoded  data.  The  combination  of  lens  LI  and  a  scene 
in  its  front  focal  plane  models  a  distant  target  such  as  the  sky.  Self- 
images  of  G1  are  formed  by  the  diffracted  light  waves  that  originated  in 
the  scene, and  the  light  intensity  falling  on  the  detector  is  modulated 
by  the  laterally  moving  G2  with  velocity  v^.  We  rotate  grating  G1  and 
G2  by  angle  cp,  keeping  their  rulings  parallel.  In  actual  practice  the 
scene  would  be  rotated  by  a  dove  prism  instead  of  rotating  the  gratings. 
The  gratings  have  a  transmission  function  g(x)  =  1  +  cos(2iTx/d).  Other 
gratings  with  more  complicated  transmission  functions  like  the  Ronchi 
ruling  can  be  used,  but  they  require  considerably  more  data  processing. 

The  scene  is  spatially  and  temporally  incoherent  so  that  the 
total  intensity  at  the  detector  is  the  sum  of  the  intensities  contributed 
by  each  point  in  the  scene.  Therefore  the  response  of  the  spectrum  to  a 
single  point  in  the  front  focal  plane  of  LI,  called  the  ’’point  spread 
response",  characterizes  the  instrument.  The  total  response  to  any 
scene  is  then  a  weighted  sum  of  point  spread  responses,  where  the  weight 
given  is  equal  to  the  intensity  for  each  point  in  the  scene.  The  compu¬ 
tation  of  the  point  spread  response  of  our  instrument  follows. 
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Point  Spread  Response 


Let  a  monochromatic  point  source  of  wavelength  X  be  located  at 
(x’,y’)  in  the  front  focal  plane  of  LI.  Therefore  a  plane  wave  with 
direction  cosines  B^”yVf  and  y  “  /l  -  ^  will  impinge 

upon  Gl.  Just  behind  G1  for  cp  =  0  we  have 

uCxjy,©"*”)  =  exp{ik(ax  +  3y)}  [1  cos(27Tx/d)  ] ;  k  =  27t/X. 

Using  a  reasoning  like  Edgar’s,  the  grating  diffracts  the  light  so  that 
three  plane  waves  propagate  beyond  the  grating  and  the  field  becomes 

u(x,y,z)  =  exp{ik(ax  +  3y  +  yz)}  + 

Y  exp{ik[(a  -  X/d)x  +  3y  +  z/1  -  (a  X/d)2  -  3^]}^“ 

Y  exp{ik[(a  +  X/d)x  +  3y  +  z/1  -  (a  +  X/d)2  -  g^]}. 

Great  simplification  results  when  the  square  root  terms  are  expanded  in 
a  power  series  and  only  the  first  three  terms  are  kept.  Since  the  instru¬ 
ment  will  be  operated  only  in  the  region  (a,B)  where  the  first  two  terms 
are  significant,  only  three  terms  are  kept,  where  the  third  term  will 
define  the  restricted  region.  Therefore, 

/I  -  (a  ±  X/d)^  -  3^  «  1  -  |<a  ±  X/d)2  _  i  ^2  _ 

^(a  ±  X/d)2  +  32]  2. 

Substituting  this  into  our  field  equation  and  collecting  like  terms 
yield 

u(x,y,z)  »  exp{ik(ax  +  gy  +  yz)  }  {l  +  exp{-(iirzX/d2) 

(1  -  3a2/2  -  32/2  -  x2/4d2)}} 

cos{(2-rr/d)  [x  -  za(l  +  a2/2  +  32/2  +  x2/2d2)]}. 
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By  restricting  the  size  of  the  source  region  the  propagating  field  is 

u(x,y,z)  ^  1  +  e  cos[27t(x  -  za)/d],  (1) 

where  the  overall  phase  factor  is  dropped  since  it  is  intensity  that  is 
detected,  and  the  approximation  is  valid  if  7rzX(3a^  +  3^  ~  X^/2d^)  <  tt/4 
and  7Tza(a^  +  3^  +  A^/d^)/d  7r/4.  The  former  condition  implies  that  the 

useful  region  (a, 3)  lies  within  an  ellipse  whose  major  axis  coincides 
with  the  3  axis  and  is  /J  times  as  large  as  the  minor  axis.  Since  the 
gratings  will  be  rotated,  the  region  of  restriction  lies  within  a  circle 
whose  radius  is  given  by  the  minor  axis,  r  =  d^/Az.  The  second  condition 
is  less  restrictive. 

The  field  behind  the  second  grating  is  given  by  the  product  of 
the  propagating  wavefield,  eq.  (1) ,and  the  grating  transmission  function 
for  G2, which  is  moving  laterally  with  velocity  v^*  Therefore 

u(x,y,z+)  {1  +  e  cos[2tt(x  -  za)/d]} 

{1  +  cos[277  (x  -  V2^t)/d]}, 

where  indicates  the  plane  just  behind  the  second  grating.  The  light 
in  this  plane  is  gathered  and  recorded  by  a  detector  whose  signal  is 
proportional  to  the  total  intensity.  Therefore  the  detected  intensity  is 
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(2) 


Ij)(z,a)  »  9/4  +  2  cos(uzX/d^)  cos[2iT(Vxt  -  za)/d] 


+  -T-  cos[4Tr(vxt  -  za)/d], 
o 


where  I  |  dxdy  is  chosen  to  be  1. 

4" 


Demodulation  of  the  Response 

The  point  spread  response  of  equ*  (2)  contains  a  dc  term  and  two 
time  modulated  terms.  Each  modulated  term  has  its  own  temporal  carrier, 
Vx/d  and  2vx/d,  and  they  can  be  filtered  either  by  a  computational  step 
or  electronically.  Therefore  the  two  terms  cos(7rzA/d^)  cos(27Tza/d)  and 
cos(47rzX/d)  can  be  separately  treated  as  point  spread  responses  of  a  new 
system  that  includes  a  filtering  step.  The  former  term  contains  both  X 
and  a  information  whereas  the  latter  has  only  a.  We  select  the  a  only 
term  for  our  Talbot  image  scanner  and  the  other  term  for  the  Talbot 
spectrometric  imager.  Before  considering  these  two  instruments  a  review 
of  the  Talbot  spectrometer  would  be  appropriate.  The  spectrometer 
selects  the  term  in  X. 


The  Talbot  Spectrometer  Review 


The  Talbot  spectrometer^  uses  the  cos (TTzX/d^)  cos(27Tza/d)  term. 
By  restricting  the  source  with  a  narrow  slit  so  that  cos(27Tza/d)  >  0.7, 
this  part  of  the  term  can  be  ignored.  Our  restriction  means  that  the 
angular  width  of  the  slit  is  a  <  d/8z.  The  detected  intensity  is 


Id(z) 


^  S(X)  cos  (iTzX/d^)  dX 


Re  S(p), 
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where  p  =  z/2d^,  Re  means  the  real  part  of  what  follows,  and  S(p)  is  the 
Fourier  transform  of  S(X).  The  spectrum  can  be  recovered  by  another 
cosine  transformation. 


cos  (27TP  A) 


dp 


j  [sa)  +  s(-A)]. 


Since  we  encode  only  the  real  part  of  the  Fourier  transform  the  recovered 
spectrum  is  even.  However  S(A)  =  0  for  A  <  0,  and  there  is  no  ambiguity 
in  the  result. 

The  resolution  of  the  spectrum  is  easily  calculated  by  letting 
S(A)  =  6(A  “  Aq  +  AA/2)+  6 (A  -  Aq  -  AA/2).  The  recorded  intensity  is 

Ij)(z)  cos[27rp(Ao  -  AA/2)]  +  cos[2'iTp(Ao  +  AA/2)] 

=  2  cos(27TpAA)  cos(47rpAo), 

where  cos(47TpAo)  is  the  carrier  and  cos(27rpAA)  is  the  envelope  of  the 

signal.  A  contrast  reduction  of  0.7  in  I]}(z)  can  be  detected  by  the  eye 

and  this  occurs  when  27rpAA  =  7r/4.  Solving  for  AA  we  have  AA  =  d^/4z. 

The  minimum  detectible  separation  6A  occurs  for  the  largest  z.  From 

chapter  II  the  largest  z  is  z  =  Nd^/2A,  where  N  is  the  number  of 

max 

grating  periods  in  the  illuminated  field.  Therefore  the  resolution 
6A/A  =  1/2N.  For  a  sodium  doublet  Aq  =  5893  A  and  AA  =  6  A  and  therefore 
a  grating  with  at  least  500  lines  can  resolve  the  doublet.  With  photo¬ 
electric  measurement  of  I(S)  and  hence  S(p)  the  resolution  will  be 
better  since  the  drop  of  the  envelope  will  be  detectible  long  before  it 
goes  down  to  0.7  of  its  maximum. 
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The  Talbot  Fourier  Analyzer 


This  instrument  uses  the  cos(47rza/d)  term  so  that  a  mask 
M(xVf>yVf)  placed  in  the  front  focal  plane  of  lens  LI  has  a  recorded 
intensity 

Id(z)  = 

where  the  interval  of  integration  is  and  a  ^  -x/f,  3  ^  -y'/f*  If 

a  real  function  g(x'/f)  is  to  be  analyzed,  then  we  make  a  function  mask 
such  that  M(xVf,y7f)  =  1  for  0  <  xVf  ^  Q,  0  <  yVf  ^  g(xVf),  and  zero 
otherwise,  where  Q  includes  the  entire  range  of  g(xVf)*  The  size  of 
the  mask  is  limited  to  the  region  where  our  simplified  point  spread 
response  is  valid.  Consequently, 


M(xVf,yVf)  cosCAttzx Vfd)  dxMy*, 


1d(z)  /  /M(xVf,y'/f)  cos(47rzx Vfd)  dxMy ’ 


dy'  cos(47Tzx Vfd)  dx* 


g(xVf)  cos(47tzx Vfd)  dx’  ==  Re  f  (v)  , 


where  v  =  2z/fd. 

Since  only  positive  values  of  g(xVf)  can  be  simulated  with  our 
mask,  functions  with  negative  values  cannot  be  directly  encoded.  A 
positive  constant  value  is  added  to  the  function  so  that  it  is  every¬ 
where  positive.  Knowing  the  result  in  the  transform  domain  of  adding 
a  constant  to  a  function,  this  bias  can  be  subtracted  from  the  final 
output . 
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The  Talbot  Image  Scanner 


Theory 

Let  S(a,3)  be  the  brightness  of  a  scene  as  a  function  of  its 
angular  coordinates  (a,  3)  and  filter  out  the  cos(47rza/d)  term  of  the 
point  spread  response*  The  recorded  intensity  will  be 


Id(z) 

where 

S(v,y) 


(3) 


is  the  two  dimensional  Fourier  transform  of  S(a,3),  v  =  z/2d  radians“^. 

The  recorded  signal  is  a  special  case  of  the  two  dimensional 

Fourier  cosine  transformation  and  explores  a  single  straight  line  in 

this  space.  The  highest  recorded  frequency  is  given  by  the  greatest 

separation  between  the  gratings  which  cannot  exceed  z  =  Nd^/2A,  where 

max  ’ 

N  is  the  number  of  grating  periods.  Another  straight  line  can  be 
explored  in  the  transform  space  by  tilting  both  G1  and  G2  by  <p.  Because 
of  this  tilt  equ.  (3)  takes  on  new  coordinates  a’,3^,v’,]j^  where  a’  = 
a  cos  <p  +  3  sin  9,  3*  =  -a  sin  cp  +  3  cos  9,  =  v  cos  9  +  y  sin  9  and 

y'  =  0.  Equation  (3)  becomes 


1d(2,9)  =  //s(a-,3’)  cos(277V*a*)  daM3*  =  Re  S(v’,0), 


where  v’  is  a  straight  line  in  the  transform  space  (v,y)  tilted  by  an 
angle  9.  This  is  sketched  in  Fig.  23. 

Thus  the  total  transform  space  is  accessible  to  the  image  scanner 
for  encoding.  The  grating  separation  z  determines  the  magnitude  of  the 
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Figure  23.  The  transform  space  for  tilted  gratings.  The  recorded  signal 
gives  us  the  value  of  the  cosine  transformation  on  the  axis  that  tilts 
with  the  gratings.  Therefore  the  entire  transform  space  can  be  sampled. 

frequency  in  the  transform  plane  that  we  are  encoding  the  tilt  of  the 

gratings  gives  us  its  orientation  so  that  the  coordinates  of  S(v,ij)  are 

given  by  V  =  v’  cos  9  and  p  =  v’  sin  9.  For  the  straight  line  sampling 

the  gratings  are  rotated  by  A9  and  then  the  gratings  are  separated  from 

zero  to  z  .  The  gratings  are  rotated  by  another  A9  and  then  separated 
max 

again.  This  is  repeated  and  in  this  manner  the  transform  domain  is 

sampled  along  straight  lines  that  look  like  a  collection  of  spikes. 

Another  way  to  sample  the  transform  space  is  by  incrementing  the 

separation  of  the  gratings  by  some  fixed  distance  Az  and  rotating  them 

through  one  revolution.  Then  they  are  separated  by  another  Az  and  once  more 

rotated.  The  sequence  is  continued  until  z  is  reached.  In  this  way 

max 

the  transform  space  is  sampled  in  regions  that  form  a  set  of  equally 
spaced  concentric  circles. 

A  third  way  to  sample  the  transform  space  is  to  separate  the 
gratings  continuously  while  rotating  them  so  that  the  transform  is  sampled 
on  a  spiral.  The  shape  of  the  spiral  depends  upon  the  relative  motions 
between  rotation  and  separation.  A  spiral  with  many  revolutions  occurs 
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if  we  rotate  many  times  as  we  separate  the  gratings.  This  situation  is 
approximately  the  same  as  the  concentric  circle  case.  A  sequence  of 
spirals  that  do  not  complete  one  revolution  can  be  made  by  separating  the 
gratings  through  the  maximum  extent  for  a  small  rotation  that  is  but  a 
fraction  of  one  revolution.  The  next  spiral  is  generated  through  separ¬ 
ating  the  gratings  again,  and  so  forth.  A  sequence  of  spirals  that  look 
like  spikes  is  obtained.  Sampling  in  the  Fourier  domain  means  that  our 
recovered  image  will  be  replicated.  If  the  sampling  is  sufficiently 
dense  then  the  total  transform  is  known  through  interpolation  in  a  way 
analogous  to  the  dimensional  sampling  theorem  in  cartesian  coordinates. 

In  this  way  the  image  can  be  recovered.  Work  by  Barakat  can  be  used 
for  the  circular  sampling  case  and  that  by  Bracewell^^  for  the  straight 
line  sampling.  Unfortunately  the  formulation  does  not  lend  itself  to 
the  Fast  Fourier  transform  for  which  the  sampling  points  should  be  on  a 
rectangular  array. 

A  less  efficient  way  of  scanning  the  transform  domain  is  to 
program  z  and  9  motions  so  that  the  scan  travels  along  lines  that  are 
parallel  to  the  v  or  y  coordinate.  This  is  accomplished  by  holding  a 
constant  value  of  say  vq  Po  depending  upon  the  scan  desired,  so  that 
Vq  =  zCt)  cos  cp(t)  and  pq  =  z(t)  sin  9(t),  where  z(t)  =  zg  +  v^t  and 
9(t)  =  9o  v^  is  the  longitudinal  velocity  of  G2,  o)  is  the  angular 

velocity  of  the  gratings  and  zq,  9q  is  the  initial  polar  coordinate  from 
which  the  scan  begins.  In  this  way  the  Fourier  domain  is  scanned  in  a 
rectangular  array  and  the  Fast  Fourier  transform  can  be  used  to  recover 
the  imag  e . 
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Special  Cases 

In  this  section  we  discuss  special  scenes  which  can  be  handled 
simply  without  the  need  for  sampling  the  transform  domain.  Consider  an 
object  that  contains  a  number  of  straight  lines  as  shown  in  Fig.  24. 


Figure  24.  Spike  scene.  This  scene  can  be  easily  deduced  by  the  instru¬ 
ment,  as  explained  in  the  text. 

The  Talbot  image  scanner  can  be  operated  in  the  following  mode  to  detect 

the  orientation  of  each  spike.  Position  G2  in  one  of  the  self-image 

planes  of  Gl.  G2  need  not  be  moved  laterally  but  can  instead  by  tilted 

slightly  with  respect  to  Gl  so  that  moire  fringes  will  be  observed  if 

there  is  a  self-image  of  Gl  at  G2.  If  this  shift  exceeds  d/4  then  the 

sum  of  the  self-images  produce  a  relatively  constant  illumination. 

Therefore  no  self-images  are  formed.  This  situation  occurs  when  the 

line  source  is  not  parallel  to  the  rulings.  By  observing  the  contrast 

of  the  moird  fringes  as  the  gratings  are  rotated  the  orientation  of  each 

line  source  can  be  determined.  From  Appendix  1,  the  source  requirement 

is  a  <  d/4z.  If  the  gratings  are  separated  by  z  =  Nd^/2A  then  the 

max 

smallest  deviation  6a  will  be  6a  <  A/2Nd.  To  convert  this  to  accuracy 
in  rotation  69,  let  L  be  the  length  of  the  source;  then  69  »  6a/L  = 
A/2NdL,  for  small  angle  approximation. 
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Another  special  object  would  be  a  uniformly  bright  scene  such  as 
a  mountain  range  at  a  far  distance  or  a  cloud  mass  of  smog.  We  use  the 
Fourier  analyzer  mode  of  the  image  scanner,  and  what  we  decode  later  on 
gives  us  the  profile  of  the  mountain  range.  If  it  is  not  uniformly 
illuminated, then  what  is  recovered  is  the  effective  profile.  This 
silhouette  example  is  mathematically  equivalent  to  the  mask  case  dis¬ 
cussed  previously  in  the  Fourier  analyzer  section. 

For  a  scene  S(a,3)  =  S^Ca)  82(3)  where  the  function  is  separable 
as  for  example  across  one  dimensional  gratings,  only  two  measurements  are 
required.  The  first  measurement  is  made  with  9=0  and  the  gratings  are 
separated.  Then  the  gratings  are  tilted  to  9  =  Tr/2  and  the  measurement 
repeated.  We  have  for  9  =  0 


82(6)  d3  I  81(a)  cos(27rva)  da 

and  for  9  =*^71/2 
r 

Si  (a)  da  I  82(3)  cos(27rv3)  d3 


C2  §1  (v) 


Cl  Re  S2(y) » 


where  Ci  =  Re  §1(0)  and  C2  =  Re  82(0).  Hence  C2  Re  Si(0)  =  C2C1  and  we 
can  use  this  to  get  Re  Si(v)/Ci  and  Re  S2(v)/C2.  The  two  dimensional 
cosine  transform  is  Re  S(v,y)  =  Re  Si(v)  Re  S2(y)*  Taking  the  inverse 
of  this  recovers  the  scene. 
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The  Talbot  Spectrometric  Imager 


The  Talbot  spectrometric  imager  uses  the  cos(7rzX/d^)  cos(27Tza/d) 
term  that  is  filtered  from  the  point  spread  response.  Therefore  for  a 
one  dimensional  scene  whose  brightness  is  given  by  S(a,X) 
the  detected  intensity  is 


Id(z) 


cos(7TzX/d^)  cos(27Tza/d) 


dadX . 


In  order  to  have  the  same  units  for  the  brightness  variable  we  substitute 
a  x/f  in  our  equation.  Expanding  the  trigonometric  function  yields 


Id(z) 


S(x,X)  {cos[27t(pX  +  vx)]  +  cos[27r(pX  -  vx)]}  dxdX 


=  Re  S(p,v)  +  Re  S(p,-v) 


(4) 


where  S(p,v)  is  the  two  dimensional  Fourier  transform  of  S(x,y),  p  = 
z/2d^,  V  =  z/fd.  Since  v  =  2pd/f,  the  complete  two  dimensional  cosine 
transform  is  not  encoded  as  a  function  of  z  but  rather  only  on  a  straight 
line  that  makes  an  angle  0  =  arc  tan  (f/2d)  with  the  v  axis,  as  shown  in 
Figure  25.  In  order  to  sample  the  whole  transform  space, the  ratio  f/2d 
must  be  varied.  This  can  be  accomplished  by  having  a  zoom  lens  arrange- 


Figure  25.  Sampling  in  the  cosine  transform  space  S(a,X).  The  sampling 
is  along  the  double  lines  which  form  an  angle  0  =  fare  tan  (f/2d). 
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merit  or  by  having  gratings  with  variable  pitch.  This  latter  arrangement 
will  be  mentioned  later. 

These  arguments  are  easily  extended  to  the  most  general  scene 
S(x,y,X).  For  this  scene  we  rotate  the  gratings  as  we  did  for  the  image 
scanner.  Therefore  variables  in  the  new  coordinate  system  become  a’  = 
a  cos  9  +  3sin  9,  3’  =  g  cos  9  -  a  sin  =  v  cos  9  y  sin  9^  y’ 

=  0;  eq .  (4)  becomes 


Ij)(z)  II  S(x,y,A)  cos{27r[pA  ±  v'(x  cos  9  +  y  sin  9)]}  dxdydX 


=  Re  S(v’,0,p)  +  Re  S(-v’,0,p), 


where  v’ ,  p  are  interdependent,  v  =  v*  cos  9,  and  y  =  y’  sin  9. 

Because  of  the  interdependence  our  scan  explores  straight  lines 
in  the  three  dimensional  transform  domain.  These  lines  go  through  the 
origin.  By  varying  9  other  such  lines  can  be  explored.  The  accessible 
lines  together  form  a  two  sided  cone  with  its  center  at  the  origin.  The 
angle  of  the  cone  is  given  by  arc  tan  (2d/f).  In  order  to  scan  the 
whole  volume  the  cone  angle  must  be  varied.  This  is  accomplished  by 
changing  f  as  with  a  zoom  lens  or  by  using  variable  pitch  gratings. 

A  variable  pitch  grating  can  be  synthesized  with  two  Ronchi 
rulings  in  contact  with  each  other  so  that  moir^  fringes  are  formed.  The 
period  p  of  the  moird  fringes  is  given  by  p  =  d/2  sin(ij;/2),  where  ij;  is 
the  angle  at  which  the  rulings  meet  and  d  is  the  period  of  the  Ronchi 
rulings.  Using  the  moir§  fringes  for  G1  and  G2  of  our  instrument  we  can 
vary  the  period  by  changing  the  angle  4^.  In  this  way  the  total  volume 
can  be  measured.  This  solution  means  a  loss  of  light  through  the  instru¬ 
ment  because  four  gratings  are  used  instead  of  two.  Therefore  only  one- 
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sixteenth  of  the  light  intensity  at  the  entrance  of  the  instrument 
reaches  the  detector  instead  of  one-fourth  with  two  gratings.  This  loss 
of  light  makes  this  a  dubious  solution.  One  alternative  is  to  use  one 
amplitude  grating  and  one  phase  grating  to  obtain  the  moire  fringes.  In 
this  way  the  light  throughput  of  this  instrument  is  unaltered. 

A  great  number  of  scanning  schemes  can  be  imagined  in  three 
dimensions.  One  would  be  a  series  of  spikes  pointing  in  all  directions. 
Another  would  be  a  three  dimensional  spiral  like  the  winding  of  a  ball 
of  string.  The  exact  scanning  scheme  has  not  been  worked  out  for  the 
Talbot  spectrometric  imager. 
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CHAPTER  V 


CONCLUSIONS 

Summary  of  Results 

A  number  of  instruments  have  been  proposed  in  this  report 
together  with  their  theoretical  analysis.  Experimental  evidence  for  the 
two  shearing  interferometers  and  the  autocollimator  has  been  presented. 
These  results  show  that  the  Talbot  effect  can  be  successfully  employed 
in  the  design  of  new  optical  instruments.  In  particular  we  mention  the 
adaptation  of  our  constant  radial  shearing  interferometer  by  a  large 
eyeglass  manufacturer  for  the  testing  of  their  new  multifocal  lenses. 

At  the  outset  we  pointed  out  the  advantages  of  the  Talbot  effect 
in  the  design  of  new  instruments.  These  advantages  were  borne  out 
during  the  course  of  the  work.  Their  low  cost  and  their  insensitivity 
to  poor  quality  gratings  is  particularly  gratifying.  Our  first  exper¬ 
iment  used  gratings  that  were  severely  scratched  without  materially 
altering  the  performance  of  the  instrument.  The  Ronchi  rulings  could  be 
purchased  for  as  little  as  $3.00  per  square  inch  (250  lines/inch)  for 

sizes  up  to  2"  x  2"^,  and  $9.00  per  sq  in.  (300  lin/in)  for  sizes  up  to 
2 

80**  X  80**  .  What  continually  amazed  visitors  about  our  demonstrations 
was  the  ease  with  which  the  setup  could  be  adjusted.  The  entire  inter¬ 
ferometer  could  be  assembled  and  adjusted  within  fifteen  minutes.  This 
is  due  to  two  factors:  (1)  The  scale  of  the  experiment  is  in  units  of 
the  self-image  interval  which  was  about  3  cm  in  our  experiments.  This 
is  a  convenient  distance  for  an  experimentalist  to  handle  since  it  does 
not  require  micropositioning  nor  an  enormously  large  laboratory. 
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(2)  The  self-image  distance  is  dependent  upon  the  grating  period.  If 
the  gratings  are  inadvertently  tilted,  then  the  effective  grating  period 
is  reduced  by  only  a  cosine  factor.  Thus  small  error  angles  of  as  much 
as  five  degrees  were  tolerable. 

The  Talbot  shearing  interferometers  measure  both  the  radial  and 
lateral  phase  derivative  of  the  wavefront  under  test  by  using  either  a 
circular  grating  or  a  Ronchi  ruling  respectively.  With  simple  binary 
filters  either  triple  or  double  shearing  Interferences  could  be  obtained. 
These  interferences  become  the  second  and  first  derivatives  respectively 
if  the  shear  introduced  is  less  than  the  width  of  the  tiniest  details 
in  the  object.  The  amount  of  shear  introduced  depends  upon  the  distance 
that  the  object  is  placed  from  the  second  grating  as  well  as  the  dif¬ 
fraction  angle.  This  affords  a  great  range  in  sensitivity.  Maximum 
contrast  in  the  interf erogram  is  obtained  by  placing  the  second  grating 
in  the  self-image  of  the  first.  Small  deviations  from  this  position 
only  affect  the  contrast  by  a  cosine  factor.  Therefore  positioning  of 
the  gratings  was  not  critical,  making  this  instrument  very  simple  to  set 
up.  Used  like  a  Schlieren  setup,  the  interferometer  can  study  the  flow 
of  hot  gases.  Beautiful  color  effects  may  be  introduced  to  give 
additional  information.  The  autocollimator  is  a  special  case  of  the 
Talbot  interferometer.  It  tests  the  curvature  of  the  wave  at  the  first 
grating  which  can  be  related  to  the  amount  of  defocussing  of  the  coll¬ 
imating  objective.  With  photoelectric  measurement  the  sensitivity  of 
the  instrument  can  be  increased  by  about  two  orders  of  magnitude. 

The  theory  of  the  image  scanner  has  been  presented,  and  special 
objects  investigated.  The  output  of  the  instrument  gives  a  single 
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point  in  the  Fourier  cosine  transform  space  of  the  scene.  By  rotating 
the  instrument  and  separating  the  two  gratings  this  point  can  be  shifted 
and  the  domain  scanned.  Several  scanning  schemes  are  discussed.  The 
Fourier  analyzer  is  a  special  case  of  this  instrument.  The  instrument 
is  not  rotated  and  the  scene  is  replaced  by  a  function  mask  in  the  front 
focal  plane  of  the  collimator.  The  preparation  of  the  function  mask  was 
described.  By  separating  the  gratings  the  cosine  transformation  of  the 
function  is  measured. 

A  theoretical  study  of  the  spectrometric  imager  is  presented. 

The  output  is  a  Fourier  cosine  transform  of  the  scene  as  a  function  of 
its  angular  coordinates  and  wavelength.  The  output  is  a  single  point  in 
the  three  dimensional  Fourier  cosine  transform  space  for  each  angular 
orientation  of  the  instrument  and  the  grating  separation.  With  two 
gratings  only,  a  surface  in  the  three  dimensional  transform  space  can  be 
scanned.  To  scan  the  complete  volume,  the  system  must  be  modified  by 
including  a  zoom  lens  or  by  using  variable  pitch  gratings.  The  variable 
pitch  grating  can  be  synthesized  by  using  the  moire  effect  between  two 
Ronchi  rulings. 


Suggestions  for  Further  Study 

It  is  the  nature  of  research  that  one  is  never  completely  done 
with  his  task.  Another  question  lurks  in  the  shadows  ready  to  replace  a 
question  recently  disposed.  In  the  following  paragraphs  are  discussions 
in  the  nature  of  unfinished  business  and  some  side  adventures. 

The  maximum  sensitivity  of  the  Talbot  instruments  is  limited  by 
the  maximum  distance  over  which  the  Talbot  images  occur.  This  distance 
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is  given  in  Chapter  II  as  Zj^  =  Bd/2X,  where  B  is  the  width  of  the 
grating,  d  is  the  grating  period  and  X  is  the  wavelength  of  the  light. 

A  Talbot  image  is  formed  by  the  interference  of  the  diffracted  orders 
with  the  zeroth  order.  These  interferences  occur  until  the  diffracted 
orders  walk  off  the  zeroth  order.  With  the  addition  of  two  mirrors 
placed  alongside  the  zeroth  order,  the  diffracted  orders  are  reflected 
back  into  the  zeroth  order  thereby  increasing  z^^  by  the  length  of  the 
mirror.  The  mirror  requirements  should  be  relatively  low  since  they  are 
used  at  a  grazing  incidence  of  the  order  of  the  diffracted  angle 
(about  lO"^  radians  for  our  experiments) .  Therefore  experiments  should 
be  conducted  to  test  this  idea.  It  would  also  be  useful  to  examine 

whether  a  number  of  smaller  mirrors  could  replace  a  single  large  mirror 

to  extend  z  „. 

MAa 

Another  interesting  experiment  suggested  by  A.  Lohmann  is  to  use 
a  prism  prior  to  the  first  grating  of  the  lateral  shearing  interferometer. 
White  light  is  used  as  a  source,  and  the  prism  diffracts  the  light  into 
its  various  colors.  Since  each  wavelength  enters  the  first  grating  at  a 
different  angle,  the  self-images  for  each  color  will  be  shifted.  This 
modified  scheme  has  the  advantage  of  tainting  the  positive  from  the 
negative  derivative  differently  since  the  setup  is  no  longer  symmetric. 

The  Fourier  analyzer,  image  scanner, and  spectrometric  imager 
should  be  constructed  and  checked  experimentally  for  feasibility. 

It  would  be  best  to  build  them  in  this  order  since  experimental  experience 
from  one  can  be  used  in  the  construction  of  the  next.  Theoretical 
studies  of  the  efficiency  for  various  scanning  schemes  on  the  latter  two 
instruments  should  be  done  so  that  an  optimum  data  collection  scheme  can 
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be  programmed.  My  guess  is  that  the  spiral  scan  is  the  most  efficient 
way  of  operating  the  instruments;  however,  recovery  of  the  scene  from 
this  type  of  scan  would  have  to  be  developed.  Schemes  for  presenting 
the  scene  as  a  function  of  its  three  parameters,  two  position  coordinates 
was  wavelength,  should  be  studied.  A  solution  would  be  to  use  a  two 
dimensional  imaging  device  such  as  a  TV  screen  and  use  color  to  present 
the  third  parameter. 

In  conclusion  we  can  say  that  our  goal  has  been  partially  accom¬ 
plished.  We  set  out  to  measure  the  distribution  of  radiation  S(a,3,l)  as 
a  function  of  its  angular  coordinates  and  wavelength.  Simple  experiments 
were  attempted  at  first  dealing  with  S(a).  From  these  studies  we  learned 
the  critical  parameters  of  the  instrument  that  we  choose.  In  particular 
we  found  that  collimation  and  good  quality  substrates  for  our  gratings 
were  important.  These  discoveries  pointed  the  way  to  the  invention  of  a 
new  autocollimator  and  two  new  shearing  interferometers.  Our  final  step 
was  to  complete  the  theory  for  an  instrument  that  would  measure  the  total 
distribution  S(a,g,X). 
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The  moir^  effect  is  very  sensitive  in  detecting  small  differences  between  two  similar  gratii^s.  These 
differences  might  be  caused  by  an  object  with  phase  gradient  placed  between  the  two  Skatings.  The  perform¬ 
ance  of  such  an  instrument  can  be  understood  in  terms  of  the  "Talbot  effect"  (also  called  . 

imaging"  or  "self-imaging").  Slight  modifications  provide  shearing  interferences  and  the  second  deriva 

tive  of  the  object. 


The  moire  effect  reveals  very  small  imperfec¬ 
tions  of  two  gratings  placed  on  top  of  each  other. 
Lord  Rayleigh  used  this  effect  for  testing  dif¬ 
fraction  gratings.  The  same  basic  idea  has  been 
used  for  many  other  purposes.  For  example,  in 
electron  microscopy  the  imperfections  in  two 
pieces  of  crystal  lattice  can  be  made  visible  in 
this  way  [1].  Essentially  the  same  moire  effect 
is  utilized  for  studying  the  shape  of  a  diffusely 
reflecting  surface,  when  the  shadow  of  a  grating 
falls  onto  that  surface.  The  grating  shadow  is 
observed  through  the  same  grating  [2].  Also  the 
shape  of  a  refracting  object  can  be  investigated 
by  means  of  the  moire  effect  [3].  The  refracting 
object  is  placed  before  or  behind  a  first  grating 
(fig.  1).  The  shadow  of  the  first  grating  will  be 
deformed  due  to  the  refractive  gradient.  The 
moire  fringes  observed  behind  the  second  grating 
placed  at  a  distance  z  from  the  object  are  lines 
of  equal  deviation  [eq.  (1)]. 


Fig.  1.  The  Talbot  fnterferometer.  Gl,  G2  gratings 
with  period  d\  ray  deflection  angle  €  in  object  O. 


e^(x,y)  =  md/z;  m  =  0,  ±1,  ±  2. . .  .(1) 

We  intend  to  extend  this  method  of  observing 
refractive  gradients  by  means  of  moire.  As  de¬ 
scribed  so  far  this  method  is  based  entirely  on 
rays.  This  point  of  view  is  not  satisfactory  since 
it  is  known  that  wave  optical  color  effects  occur 
when  white  light  passes  through  two  gratings  at 
a  finite  distance.  This  happens  with  uncollimated 
[4]  and  with  collimated  light  [5].  Such  color  ef¬ 
fects  may  well  prove  to  be  useful.  A  more  im¬ 
portant  objection  against  the  ray-optical  point  of 
i  view  is  the  inability  to  explain  the  fundamental 
limitation  of  this  method.  Based  on  eq.  (1),  one 
j  might  suspect  that  an  arbitrarily  small  deflec¬ 
tion  angle  can  be  detected  if  only  the  distance  2 
of  the  two  gratings  is  sufficiently  large.  Since 
such  unlimited  detectibility  is  never  attainable 


Fig.  2.  Overlapping  diffraction  orders  behind  a  grating 
of  finite  width  B  =  Nd. 
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(b) 

Fig.  3.  Talbot  interferences  with  a  candle  flame  as  object  between  the  two  gratings,  a.  Gratings  parallel,  b.  One 

grating  slightly  rotated. 
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one  might  try  to  patch  on  the  following  wave-op¬ 
tical  argument. 

The  light  passing  through  a  grating  slit  of 
width  d/2  does  not  behave  like  a  ray  anymore 
after  it  has  travelled  over  a  distance  where 
the  diffraction  spread  Zj^X/(d/2)  equals  the  slit 
width  d/2  [eq.  (2)]. 

z^X/{d/2)  =  d/2  ;  z^  =  d^/AX  ,  (2) 

Based  on  this  criterion  one  might  expect  to 
be  the  smallest  detectable  deflection  [eq.  (3)]. 

=  d/2z^  =  2X/d  .  (3) 

With  X  =  0.5  X  10"^  mm  and  =  10"^  this  leads 
to  Cp  =  which  is  only  moderately  good. 

Fortunately  the  method  is  capable  of  detecting 
angles  much  smaller  than  Cp  [eq.  (3)].  This  is 
due  to  the  Talbot  effect  [6],  which  is  also  known 
as  "Fourier  imaging"  [7]  and  "self-imaging"  [8]. 
Talbot  discovered  about  135  years  ago  that 
images  can  be  formed  without  any  lenses  or 
mirrors  if  the  object  is  a  grating  which  is  il¬ 
luminated  in  collimated  monochromatic  light. 
These  "Talbot  images"  occur  at  distances  2d^/X, 
4c?2/X,  etc.  behind  the  grating.  When  a  second 
grating  at  a  slight  angular  rotation  to  G1  is 
placed  into  the  plane  of  a  Talbot  image,  moire 
fringes  of  high  contrast  are  observed.  These 
moire  fringes  will  be  deformed  if  a  refractive 
phase  object  is  placed  for  example  close  to  the 
first  grating.  As  before  these  moire  fringes  in¬ 
dicate  lines  of  equal  deflection  by  the  object 
[eq.  (1)]. 

The  largest  possible  distance  Zj,  (or  maybe 
2>2:  j)  depends  on  the  finite  width  B  =  Nd  of  the 
first  grating.  At  the  first  grating  diffraction 
orders  have  moved  to  both  sides  by  half  of  the 
grating  width  B  [fig.  2;  eq.  (4)]. 

Zrj.x/d  =  B/2  =  Nd/2  ;  zt  =  Nd^/2\  .  (4) 

When  the  two  gratings  are  separated  by  the 
smallest  detectable  deflection  Crp  is  now  smaller 
by  a  factor  2^^,  where  N  is  the  number  of  periods 
in  the  first  grating  [eq.  (5)]. 

€t  =  d/2zrY  -  X/Nd  =  X/B  =  €p/2iV^  .  (5) 

We  have  performed  some  experiments  with 
this  "Talbot  interferometer”  as  shown  in  fig.  1. 

In  fig.  3a  the  object  is  a  flame,  and  again  in 


fig.  3b,  but  now  with  the  second  grating  slightly 
rotated  around  the  optical  axis.  In  some  addi¬ 
tional  experiments  we  have  placed  a  telecentric 
lens  system  behind  the  Talbot  interferometer. 
When  introducing  a  spatial  filter  into  this  tele- 
centric  system  and  rejecting  everything  but  the 
zeroth  diffraction  order  from  the  two-grating 
system  one  observes  in  essence  the  second 
derivative  d^u/dx"^  of  the  object  w(Ar,y).  When 
shifting  the  spatial  filter  to  one  of  the  first 
grating  diffraction  orders  the  image  represents 
shearing  interferences  u{x  +Xz/d^  y)  -  uipc^  y). 

In  a  white  light  one  obtains  beautiful  color 
fringes  which  are  unlike  ordinary  interference 
fringes.  For  example  for  a  specific  grating 
distance  the  image  contains  many  blue -orange 
fringes.  These  color  fringes  are  of  such  high 
contrast  that  E.  Lau  found  a  similar  setup  useful 
in  designing  tablecloths  [9].  Our  method  is  also 
applicable  for  the  detection  of  small  differences 
between  two  quite  irregular  but  similar  objects. 
One  begins  by  recording  photographically  the 
fringes  from  the  first  object.  After  development 
the  photograph  is  placed  where  the  fringes  had 
been  observed.  The  first  object  is  now  replaced 
by  the  second  object,  which  might  actually  be 
the  first  object  but  somewhat  deformed.  The 
moire  fringes  between  subsequently  produced 
Talbot  interference  fringes  will  reveal  small 
differences  between  the  two  objects.  The  quanti¬ 
tative  evaluation  is  similar  to  that  for  Lau^s 
dupligram  method  [10]  and  for  life-fringe  holo¬ 
graphic  interferometry. 
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The  linear  gratings  in  our  Talbot  interferometer  are  replaced  by  circular  gratings.  This  new 
radial  shearing  interferometer  is  sensitive  to  the  radial  gradient  of  the  object. 


Talbot  [1]  discovered  in  1836  that  images  of  a  photographic  reduction.  The  interference  image 
grating  can  be  formed  without  any  lens  or  mirror,  of  a  bifocal  glass  is  shown  in  fig.  2.  The  inter- 
When  illuminated  by  a  monochromatic  plane  wave  ference  fringes  degenerate  into  spikes  (actually 
the  grating  will  be  imaged  at  distances  stretched  hyperbolas)  when  the  two  circular 

etc.  (d  =  grating  constant).  Recently  this  gratings  are  not  laterally  aligned.  These  spikes 
Talbot  effect  has  been  used  for  interferometry  [2],  assume  characteristic  shapes  when  a  phase  ob- 
A  second  identical  grating  is  placed  into  one  of  ject  such  as  a  candle  is  introduced  (fig. 3). 
the  "Talbot-planes"  of  the  first  grating.  The  There  are  some  remaining  questions,  which 

phase  object  under  investigation  is  situated  be-  we  intend  to  answer  in  a  more  elaborate  paper 

tween  the  two  gratings.  The  image  consists  es-  in  "Applied  Optics".  For  example,  when  this 

sentially  of  two  shifted  object  wavefronts  or  of  circular  Talbot  interferometer  is  to  be  used  for 
the  first  derivative  du(x,y)/dx  of  the  object  quantitative  investigations  one  has  to  study  how 

u(x,y).  well  the  circular  grating  really  produces  "  self- 

For  some  objects  the  radial  derivative  images"  [3].  In  connection  with  a  spatial  filter 


du(r,cp)/dr  is  more  indicative  than  the  linear 
derivative  du/dx.  We  obtain  this  radial  gradient 
by  replacing  the  linear  gratings  by  radial  gratings 


I  in  plane  F  of  fig.  1  this  interferometer  can  also 
I  produce  the  second  derivative  d'^u/dr^.  Or  it 


as  shown  in  fig.  1.  Our  radial  gratings  are  from 


a  moire  kit*  with  a  radius  step  of  0.25  mm  after  *  Product  of  Edmund  Scientific,  Barrington,  N. J. 
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Fig.  2.  Interference  image  of  a  bifocal  eyeglass.  The  gratings  G1  and  G2  are  centered. 


Fig.  3.  Interference  image  of  a  candle.  The  two  gratings  are  not  centered. 
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There  are  two  common  methods  of  determining  the  degree  of 
collimation;  autocollimation  and  shearing  interferometry. 
The  simplest  method  is  the  autocollimation  technique,  but  it 
only  indicates  collimation  by  comparing  the  size  of  the  source 
with  its  image.  On  the  other  hand  the  spacing  of  the  fringes  in 
shearing  interferometry  is  a  direct  measure  of  the  degree  of 
collimation.  1  Recently,  Langenbeck  derived  a  method  that  used 
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two  tiny  corner-cube  reflectors  to  sample  the  beam.  In  this 
fashion  he  was  able  to  translate  the  measurement  from  fringe 
spacing  to  fringe  rotation,  with  the  result  of  increased  sensitivity.^ 
The  method  presented  here  is  an  application  of  the  Talbot  inter¬ 
ferometer^  and  has  the  same  limitations  as  the  other  interfero¬ 
metric  techniejues.  The  main  advantages  of  this  technique  are 
the  inexpensive  components  (two  Ronchi  rulings)  and  the  relative 
insensitivity  to  (ornponent  alignment. 

As  Talbot  observed  ini  (So7,  an  image  of  a  grating  appears  at 
integral  multiples  of  the  distance  2d^/\  (where  d  is  the  grating 
period)  when  the  grating  is  illuminated  by  a  plane  monochromatic 
wave.  A  Ronchi  ruling  would  therefore  be  imaged  at  2dV^) 
4d^/\  etc.  If  we  place  another  identical  Ronchi  ruling  in  one 
of  these  self-image  planes  (see  Fig.  1 )  moir6  fringes  will  be  formed, 
as  observed  by  J.  Burch.  ^  The  fringe  spacing  P  can  be  calcu¬ 
lated  by  using  the  vector  diagrams  described  by  Rogers  Ho  be 

P  =  d/2sin(0/2),  (1) 

where  0  is  the  angle  with  which  the  two  gratings  meet.  As  the 
grating  G2  is  rotated  around  the  optical  axis,  the  fringe  spacing 
increases  until  we  have  uniform  brightness  (or  darkness)  when 
grating  bars  of  Gl  and  G2  are  parallel  to  each  other. 

The  plane  wave  that  illuminates  Gl  of  the  interferometer  is 
obtained  by  placing  a  point  source  at  the  focus  of  a  lens  as  shown 
in  Fig.  1.  When  this  lens  is  defocused,  the  plane  monochromatic 
wave  becomes  spherical.  Cowley  and  Moodie®  have  shown  that 
the  positions  of  the  self-image  planes  of  a  grating  illuminated 
by  a  spherical  wave  are  given  by 

l/Zj  +  l/Zfv  -  1/2/Zr,  (2) 

where  «/  is  a  positive  integer,  1/Zw  is  the  curvature  of  the  in¬ 
cident  wave  at  Gl,  Zj  is  the  location  of  the  self-image  as  measured 
from  Gl,  and  Zt  is  the  distance  2d^/\.  The  self-image  will  be 
magnified  by  a  fa(!tor 

A/  =  I  +  ZjjZw.  (3) 

This  image,  when  superimposed  with  the  second  grating,  again 
produces  moire  fringes.  Now  the  grating  G2  is  rotated  around 
the  optical  axis  until  the  giating  bars  of  Gl  and  G2  are  parallel 
to  each  other.  In  this  case,  the  Rogers  vector  diagram  gives  a 
fringe  spacing  P  —  did2/\di  —  d-ajj  where  di  is  the  period  of  the 
self-image,  da  the  period  of  the  second  grating,  and  di  =  Md. 
The  fringe  spacing  in  teinis  of  the  grating  period  d  is  therefore 
P  -  Md/\M  —  l|.  Thus  the  fringe  spacing  is  a  quantitative 
measure  of  the  degree  of  collimation,  since  |M  —  1 1  is  propor¬ 
tional  to  the  curvature  l/Zw, 

Now  we  want  to  discuss  the  accuracy  obtainable,  that  is, 
the  smallest  detectable  deviation  from  perfect  collimation.  If 
the  collimation  is  perfect,  the  magnification  M  of  the  Talbot 
image  is  1,  and  hence  the  moire  period  P  is  infinite.  Thus  our 
task  is  to  find  how  small  the  period  P  must  become  in  order  to  be 
detectable.  Expenence  has  shown  that  the  presence  of  a  moir6 
effect  is  clearly  visible  if  at  least  one-half  of  the  moir^  period  P 
falls  within  the  limited  field  of  observation  with  width  B.  Thus 
one  fringe  will  be  detected  whenever  the  number  of  lines  between 
grating  G2  and  the  image  of  Gl  differ  by  one-half.  With 
divergent  illumination  we  count  N  lines  of  G2  in  the  observation 
field  of  width  B  and  A'  —  §  lines  of  the  image  of  Gl.  Thus 
the  magnification  M  required  for  fringe  detection  is  bounded  by 

M  >  N/{N  -  I)  «  1  +  1/2A,  (4) 

and  in  convergent  illumination, 

M  <  N/{N  +  i)  «  1  -  1/2A.  (5) 


Fig.  1.  Talbot  interferometer  setup  for  observing  fringes  due  to 
defocusing  of  collimator  objective  {Zo  5*^  /).  Fringes  are  observed 


Fig.  2.  The  self-imaging  effect  occurs  within  overlapping  orders. 
Test  done  at  distance  where  the  two  first-order  diffraction  spots 
begin  to  separate. 


M  -  (1  -  2JZt/ZwY^  «  1  +  2JZtIZw^  (6) 

If  M  falls  outside  the  bounds  of  Eqs.  (4)  and  (5),  no  fringes  will 
appear,  limiting  the  sensitivity  of  this  test.  Solving  Eq.  (6) 
with  the  bounds  of  Eqs.  (4)  and  (5),  we  find  that  the  minimum 
detectable  field  curvature  l/Zw  is  |l/^Tr|  >  (4Art/ZT)”^,  where 
the  absolute  signs  obviate  the  necessity  of  assigning  a  sign 
convention  for  field  curvature.  By  substituting  Eq.  (2)  into 
the  above  equation  we  obtain,  ior  N  ^  1, 

\l/Zw\  >  {2NZj)~K  (7) 

From  Eq.  (7)  it  follows  that  a  large  value  of  Zj  is  desired  in 
order  to  detect  small  collimation  defects  l/Zw-  The  farthest 
distance  to  which  we  can  go  is  limited  by  the  walkoff  of  the  first 
order  grating  diffraction.  As  indicated  in  Fig.  2,  the  longest 
distance  at  which  there  is  still  a  connected  interference  field  of 
width  D  is  at  Zj  =  Dd/2\.  Perhaps  one  could  perform  collima¬ 
tion  tests  up  to  about  twice  this  distance,  but  we  will  restrict 
ourselves  to  Zj  <  Dd/2\.  Inserting  this  limit  into  Eq.  (7), 
together  with  N  —  D/d^  we  get 

\l/Zw\  >  2\/DK  (8) 

The  curvature  of  the  wave  incident  on  Gl  is  related  to  the 
defocusing  of  the  collimating  objective  by 

\l/Zw\  =  1/  -  Z^\l(Za\S  -  ^o|  +  ^/)  «  1/  -  -Zol/^o/  (9) 

where  Zq  is  the  distance  of  the  collimating  objective  from  the 
point  source,  Zq  is  the  distance  from  objective  to  Gl,  and  / 
is  focal  length  of  the  objective  lens.  The  approximation  is 
permissible  since  the  distance  Zg  “  Zq\  in  the  test. 

Call  5  =  1/  —  Zo|  the  focusing  error,  and  substitute  Eq.  (8)  into 
Eq.  (9).  We  have 


By  combining  Eqs.  (2)  and  (3)  we  have 


5  >  2HS/D)\ 


(10) 


where  we  have  used  the  fact  Zo  ~  /.  Thus  the  focusing  error 
is  simply  related  to  the  wavelength,  and  the  square  of  the  focal 
length  to  beam  diameter  ratio. 

The  focusing  error  can  be  reduced  by  simply  rotating  <ine 
gi'ating  with  respect  to  the  other,  thereby  producing  moir4  fringes 
in  the  field  of  observation.  In  this  fashion  we  change  our  de¬ 
tection  scheme  from  the  condition  of  uniform  brightness  to  a 
system  of  rotating  fringes.  It  is  known  that  fringe  detection 
is  more  accurate  than  use  of  the  uniform  brightness  condition. 
We  rotate  both  gratings  in  opposite  directions  by  angles  6/2  with 
the  y  axis.  If  both  gratings  have  the  same  period,  moir6  fringes 
will  appear  in  the  field  of  observation  parallel  to  the  x  axis  with 
spacings  given  by  Eq.  (1).  When  the  periods  are  unequal, 
the  fringe  spacing  is  given  by 

P  =  dido/idi^  +  ^2^  —  2did2  cos0)^  —  did2/\di  —  ^2!, 

where  the  approximation  is  for  small  angular  rotations  of  0,  and 
dijd2  are  the  grating  periods  of  G1  image  and  G2,  respectively. 
It  can  be  shown  by  using  the  vector  arguments  of  Rogers  that 
these  fringes  of  unequal  grating  periods  are  rotated  through  an 
angle  with  respect  to  the  x  axis,  and  the  direction  of  this  rotation 
depends  on  whether  the  illumination  is  converging  or  diverging. 
Thus  the  experimenter  knows  in  which  direction  to  move  the 
objective  lens.  This  rotation  is  related  to  the  magnification  M  by 
M  —  (cot^/2  -f  tan<^)/(cot^/2  —  tan^).  Substituting  this 
equation  into  Eq.  (3)  we  obtain 

l/Zw  —  {l/Zj)[2  ia,n4>/{coW/2  —  tan^)].  (11) 


minus  first  zeroth  plus  first  order 


Fig.  3.  Moire  fringes  which  are  rotated  through  angle  <i>  because 
of  defocusing  of  collimator  objective. 


Our  ability  to  measure  <f>  limits  the  sensitivity  of  this  test.  If 
we  say  the  minimum  detectable  fringe  angle  0  is  about  one-half 
of  a  fringe  (see  Fig.  3),  then  Eq.  (11)  reduces  to  \  l/Zw\  >  2X/D^. 

This  lower  bound  is  identifial  to  that  of  the  uniform  brightness 
(jondition  [Eq.  (8)].  In  the  laboratory  one-quarter  of  a  fringe 
was  easily  discernible,  producing  results  which  were  a  factor  of  2 
better  than  predicted.  If  we  want  to  achieve  even  higher  ac¬ 
curacy,  we  may  use  photodetectors  instead  of  the  eye.  Assuming 
1%  brightness  accuracy,  we  can  detect  fringes  one  hundred 
times  larger  than  the  field  of  view.  Hence  the  photoelectrical 
approach  would  improve  the  detectibility  of  collimation  errors 
by  a  factor  of  about  oO  over  that  which  was  predicted. 

The  visual  experiments  have  confirmed  our  theoretical  per¬ 
formance  predictions.  An  objective  lens  of  20  cm  focal  length 
was  used  to  produce  a  5-cm  beam.  The  light  from  a  helium- 
neon  laser  was  focused  with  a  3.9-mm  focal  length  lens  onto  a  6-^* 
pinhole.  The  position  accuracy  was  0.001  cm,  whereas  the 
calculated  accuracy  was  0.002  cm.  This  improvement  was  due 
to  the  fact  that  one-quarter  fringe  was  used  with  the  angular 
measurement  technique. 

Last,  one  can  improve  the  accuracy  of  the  test  by  increasing 
the  distance  Zj.  This  distance  was  limited  by  the  walkoff  as 
shown  in  Fig.  2.  By  the  addition  of  one  or  two  mirrors  placed 
where  the  borderline  rays  of  the  zero-order  diffraction  propagate 
(Fig.  3),  the  diffracted  first-order  beam  can  be  reflected  back 
into  the  zero-order  beam  to  create  the  shearing  interferences. 
This  should  increase  Zj  and  the  accuracy  of  the  test  manyfold. 
However,  this  has  not  been  verified  experimentally. 

The  author  is  indebted  to  Adolf  Lohmann  for  his  support  and 
interest  in  this  work  and  for  his  helpful  suggestions  in  the  prepara¬ 
tion  of  this  Letter.  This  work  was  performed  under  Air  Force 
contract  AF-F19628-69-C-0268. 
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Appendix  1 


Source  Size 

A  point  source  will  give  self-images  of  the  grating  in  a  region 
where  the  diffracting  orders  interfere.  As  the  size  of  the  source 
increases,  the  distance  over  which  self— imaging  occurs  decreases  until 
images  are  no  longer  observed.  The  spatial  coherence  (size  of  the 
source)  is  derived  with  the  aid  of  Figure  27. 


Figure  26.  Diagram  showing  the  angular  size  of  the  source.  Source 
points  off-axis  produce  shifted  Talbot  images.  Summing  the  contribu¬ 
tions  from  all  points  on  the  extended  source  yields  a  self-image  of  G1 
only  if  we  do  not  allow  the  image  to  shift  by  more  than  say  a  quarter 
of  a  period,  a. 

An  off-axis  point  source  will  produce  self-images  that  are 
centered  along  the  principal  ray,  that  forms  an  angle  ot  with  the  optical 
axis,  and  at  distances  that  are  a  multiple  of  2a^/X,  where  a  is  the 
g];^3^ring  period  and  X  is  the  wavelength  of  the  light  source.  The  shift 
s  in  the  image  is  given  by  otz.  We  have  s  =  otz  =  2Ma^/X,  where  M  is  an 

integer.  Each  off-axis  point  in  an  extended  source  contributes  a 
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shifted  image  to  the  observer.  If  the  self-image  of  the  grating  is  to 
be  observed,  the  contributing  image  from  the  furthest  off-axis  point 
must  not  be  shifted  by  more  than  a  quarter  of  a  period,  i.e.  s  =  2aMa^/X 
<  a/4,  or  a  <  A/8Ma. 

The  source  size  is  related  to  a  by  a  =  Ax/2f,  where  Ax  is  the 
size  of  the  source  and  f  is  the  focal  length  of  the  collimating  lens. 
Therefore  for  a  self-image  to  occur  at  some  multiple  M  of  the  self-image 
distance  2a^/A,  the  size  of  the  source  must  be  limited  by  Ax  <  Af/4Ma. 


Appendix  2 

Approximations  Leading  to  the  Second  Derivative 


Our  expression  for  spatially  filtered  fringes  contains  a  summa¬ 
tion  of  derivatives  that  resulted  from  a  Taylor  series  expansion.  In 
this  appendix  we  show  that  the  lowest  order  term  dominates  the  series, 

i.e. 

00 

(x,y)  A^/p!  «  Uq^^(x,y)  A^/2 

p=2 


if 


|uQ^\x,y)|  A^/p!  <  n[up^(x,y)] 


max 

where  uo^\x,y)  =  3^  uo(x,y)/3x^,  [uo^^(x,y)]  is  the  maximum  value 

max 

(2) 

that  the  function  uq  (x,y)  obtains,  A  =  Z2A/d  is  the  lateral  shift  and 
n  is  a  factor  which  controls  the  degree  to  which  the  approximation  is 
valid  and  is  set  to  1/10  for  the  remaining  paragraphs.  Rearranging  this 
equation  we  get 


[up^x.y)]^^^  >  |u$P\x,y)|  aP~^  20/p!. 


(1) 
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By  the  theorem  of  Berstein,  the  derivatives  of  a  bandlimited 
function  UQ(x,y)  are  bounded  by  the  maximum  value  that  the  function 
attains,  i.e. 

luo^\x,y)l  <  (irAv)^  [uq (x,y) 

where  Av  is  the  bandwidth  of  the  function  Uo(x,y)  in  the  Fourier  domain 
and  run(x,y)]  is  the  maximum  value  of  this  function.  This  theorem 
implies  that 

or 

[up^(x,y)]  >  lu^P\x,y)l  (ttAv)^ 

maA 

Eq.  (1)  is  the  condition  for  (x,y)  under  which  our  approximation 

for  the  second  derivative  is  valid.  This  maximum  is  bounded  by  the 

lateral  shift  A  which  is  under  the  control  of  the  experimenter.  The 

bound  of  eq.  (1)  must  be  tighter  (i.e.  lower)  than  the  above  equation. 

Therefore  A^  ^  20/p!  s  (ttAv)  or  A  <  (p!/20)  /ttAv. 

The  factorial  dominates  the  root  process  as  can  be  shown  by 

using  Stirling's  approximation  for  factorial 

1  1 

(p!/20)^“^  ~  (/2Tr/20)^~^  (p/e)  for  large  p. 

Therefore  the  smallest  p  (p  =  4)  is  the  tightest  bound  for  A. 

A  <  i/^/2Q)  (4/e)/TrAv  »  6x/6 

where  6x  =  1/Av.  This  equation  states  that  the  shear  must  be  less  than 
about  one-sixth  of  the  width  of  the  object  details  in  order  that  the 
second  derivative  dominate  the  Taylor  series.  This  estimate  is  a  worst 
case  consideration.  It  is  common  practice  in  differential  shearing 
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interferometry  to  use  a  shift  of  a  half  resolution  unit  6x/2  instead  of 
6x/6  without  noticeable  detrimental  effects.  The  reason  for  using  some 
what  larger  shifts  is  improved  detectability  of  soft  slopes  in  the 
object. 


Appendix  3 

Theoretical  Consideration  for  First  Order  Filtering 


By  blocking  one*half  of  our  grating  with  an  opaque  screen  our 
grating  transmission  function  G(r)  is  modified  as  follows: 


G(r)  =/ 


^  27rimr/a  ,  *  ^ 

Cjj^e  ,  4a<r<A,  0<cp<Tr 


otherwise. 


This  in  turn  necessitates  a  reappraisal  of  our  theoretical  calculations 
Fortunately  as  we  shall  soon  see  our  modified  setup  differs  from  the 
theory  by  only  a  multiplicative  constant  that  accounts  for  the  reduced 
intensity. 

With  our  new  limits  the  field  beyond  G1  is  given  by 

•oo 


u(r,cp,z)  ~ 


u(r',zi+) 


'0  JO 


exp{iTr[r^  +  -  2rr’  cos(<p  -  cp’)]/Az}  r'dr'(p’» 


where  we  recognize  that 


‘27r 


-ixcos6  1 

e  de  =  — 


-ixcos0 
e  de . 


'0  *^0 
Therefore  the  only  difference  in  the  computation  of  this  equation  from 
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that  of  the  unblocked  grating  equation  is  the  factor  2.  Thus  we  copy  the 
results  obtained  previously  and  divide  by  2  to  obtain 


1  _  /r-  7-  7-“ - 7 — 7  iTrXzi  (m/a)^  27riinr/a 

u(r,z)  T  (mAzi/ra)  e  e 

(m) 

The  remaining  integrations  are  the  same,  so  that  behind  the  final 
grating  we  have 


u(r,cp,Z2+)  ~ 


iTrm  27ri(m+n)r/a 
e  e 


(m)  (n) 


/l  +  mMa/r  uo(r  mAz2/a,cp). 


This  result  differs  from  our  theoretical  result  for  full  circular 
gratings  by  only  a  constant  factor.  Thus  blocking  one-half  of  the 
circular  grating  only  reduces  the  overall  intensity  while  allowing  us  to 
separate  out  the  first  diffracted  orders. 


Appendix  4 

Tolerances  for  Spatial  Filtering  a 
Test  Object  with  Multiple  Focal  Lengths 

Consider  a  test  object  with  focal  lengths  f^  and  f2  corresponding 
to  test  object  areas  0^  and  O2  respectively.  The  diameter  of  0^  is  h^ 
and  of  O2  is  h2.  The  effective  focal  lengths  are  given  by 

f  =  fifL/(fl  +  fL  -  D),  f"  =  f2fL/(f2  +  fL  -  D), 
where  is  the  focal  length  of  lens  L2  and  D  is  the  distance  that  L2  is 
from  the  test  object. 

The  diffracted  orders  therefore  fall  at  distances  f  and  f"  from 
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Figure  27(a).  The  tail  end  of  the  Talbot  interferometric  setup  showing 
how  a  bifocal  lens  causes  the  zeroth  diffracted  order  to 
meet  in  two  different  planes  f’  and  f**. 


Figure  27(b).  An  enlargement  of  the  filtering  region  that  lies  between 

f’  and  f*'  in  (a).  The  plane  z  =  F  is  the  optimum  location 
for  separating  out  either  the  zeroth  or  first  diffracted 
orders,  y  +  2x  is  the  distance  between  the  zeroth  and 
first  orders. 


Figure  27(c).  An  expansion  of  the  triangular  region  in  (b)  which  aids  in 
the  computation  of  x. 
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the  principal  plane  of  the  combined  system  0  and  L2,  The  angles  at 
which  the  diffracted  orders  meet  are  a*  and  a”*  These  quantities  are 
marked  in  Figure  27(a).  An  expanded  view  of  the  area  in  which  the 
diffracted  orders  meet  is  drawn  in  Figure  27(b).  We  see  that  filtering 
can  be  accomplished  in  plane  F  as  long  as  the  diffracted  orders  do  not 
overlap.  This  is  equivalent  to  saying  that  the  distance  marked  y  must 
be  greater  than  zero.  In  Figure  27(b)  we  have  drawn  the  figure  as  though 
the  first  diffracted  order  for  0]^  lay  at  the  same  distance  from  the 
optical  axis  as  O2  (i.e.,  Af'/a^  Af"/a).  This  assumption  is  a  conser¬ 
vative  one,  since  the  distance  y  is  slightly  larger  if  we  do  not  make 
this  assumption.  To  summarize,  filtering  is  possible  at  plane  z  =  F  if 
y  =  (AfVa)  -  2x  >  0,  or  A/a  >  2x/f\ 

The  distance  x  is  computed  from  a  trigonometric  consideration  of 
the  oblique  triangle  containing  x  and  a  side  of  length  Af.  This  triangle 
is  enlarged  in  Figure  2  7(c) .  From  the  law  of  sines  we  have  a/Af  = 
sin(a*V2)/sin(a V2  +  a*V2),  and  the  right  triangle  gives  us  x  = 
a  sin(aV2),  where  sin(a'/2)  «  tan(aV2)  =  hi/2f’  and  sin(a*72)  ^ 
tan(ot*V2)  =  h2/2f".  Therefore  x  Af  hih2/2(hif"  +  h2f’).  Our  bound  on 
the  diffraction  angle  is  therefore  A/a  >  Af  hih2/ (f  ^ )  ^  (h]^  +  h2),  where 
we  have  assumed  that  f'  ^  f**.  Some  simplification  will  result  if  we  let 


1  -  i 

fi  ft. 


-2 - iju. 

fL  ^  f2  fr  ^  f2  fl  ^  ’ 


where  M  is  the  lateral  magnification.  Hence 


X/a  >  ( 


and  we  conclude  that  the  focal  length  of  lens  L2  is  immaterial  for  these 


considerations . 
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